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ABSTRACT 

Using VLT/FLAMES optical integral field unit observations, we present a detailed study of 

r ~~^ UM 448, a nearby Blue Compact Galaxy (BCG) previously reported to have an anomalously 

high N/O abundance ratio. NTT/SuSI2 images reveal a morphology suggestive of a merger of 

\^ two systems of contrasting colour, whilst our Ha emission maps resolve UM 448 into three 

• separate regions that do not coincide with the stellar continuum peaks. UM 448 exhibits com- 

'-^ plex emission line profiles, with most lines consisting of a narrow (FWHM<100 kms -1 ), 

central component, an underlying broad component (FWHM^ 150-300 kms -1 ) and a third, 
narrow blue-shifted component. Radial velocity maps of all three components show signs of 
S-j solid body rotation across UM 448, with a projected rotation axis that correlates with the con- 

tinuum morphology of the galaxy. A spatially-resolved, chemodynamical analysis based on 
the [O in] AA4363, 4959, [N n] A6584, [S n]AA6716, 6731 and [Ne in] A3868 line maps, is 
presented. Whilst the eastern tail of UM 448 has electron temperatures (T ) that are typical 

y— \ of BCGs, we find a region within the main body of the galaxy where the narrow and broad 

J> [O III] A4363 line components trace temperatures differing by 5000 K and oxygen abundances 

t-H differing by 0.4 dex. We measure spatially resolved and integrated ionic and elemental abun- 

dances for O, N, S and Ne throughout UM 448, and find they do not agree, possibly due the 
flux- weighting of T c from the integrated spectrum. This has significant implications for abun- 
dances derived from long-slit and integrated spectra of star-forming galaxies in the nearby and 

Q\ distant universe. A region of enhanced N/O ratio is indeed found, extended over a ^0.6 kpc 2 

region within the main body of the galaxy. Contrary to previous studies, however, we do not 

04 find evidence for a large Wolf-Rayet population, and conclude that WR stars alone cannot be 

responsible for producing the observed N/O excess. Instead, the location and disturbed mor- 
phology of the N-enriched region suggests that interaction-induced inflow of metal-poor gas 

• may be responsible. 

x 
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1 INTRODUCTION 

Observations of Blue Compact Galaxies (BCGs) in the nearby Uni- 
verse provide a means for studying chemical evolution and star 
formation (SF) processes in chemically un-evolved environments. 
Their low metallicities make them attractive as analogues to the 
young 'building-block' galaxies thought to exist in the high-z pri- 
mordial universe (Searle & Sargent 1972). BCGs are thought to 
have experienced only very low-level star formation in the past, 
but are currently undergoing bursts of SF in relatively pristine 
environments, ranging from 0.02-0.5 solar metallicity (Kunth & 
Ostlin 2000). Thus understanding the chemical evolution within 

* E-mail: bjames@stsci.edu 



BCGs can impact upon our understanding of primordial galaxies 
and galaxy evolution in general. 

The distribution of nitrogen as a function of metallicity in SF 
galaxies has often been a cause for debate, especially at interme- 
diate metallicities (7.6< 12+log(0/H) <8.3). Within this metallic- 
ity range the N/O abundance ratio indicates that nitrogen behaves 
neither as a primary element (i.e. N/O independent of O/H) nor 
as a secondary element (i.e. N/O proportional to O/H); instead a 
large scatter is observed. In particular, several BCGs exist whose 
N/O abundance ratios are up to three times those expected for their 
O/H values, as highlighted by e.g. Pustilnik et al. (2004). Several of 
these galaxies have been the focus of our spectroscopic studies us- 
ing integral field unit (IFU) observations (James et al. 2009, 2010), 
aimed at understanding the source of such abundance peculiarities. 
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James et al. (2009) performed a detailed abundance analysis 
of the N/O-outlier BCG, Mrk 996, whose X-ray properties have 
been studied by Georgakakis et al. (2011). Previous long-slit ob- 
servations of this galaxy (Thuan et al. 1996; Pustilnik et al. 2004), 
had shown it to have N/O ratio values at least ~4 times higher 
than those of other BCGs of similar metallicity. Using IFU obser- 
vations, James et al. (2009) performed a spatially resolved multi- 
component emission line analysis of Mrk 996, which revealed a 
twenty-fold N/O and N/H enrichment factor in the broad compo- 
nent (400-500 km/s full width at half maximum - FWHM) versus 
the narrow component line emission. The narrow component gas 
on the other hand has a N/O ratio which is typical for the galaxy's 
metallicity. In this respect, Mrk 996 is similar to NGC 5253, a star- 
burst galaxy for which an enhanced N/O ratio has also been iso- 
lated to broad component gas (Lopez-Sanchez & Esteban 2010). 
In Mrk 996, the spatial morphology of the broad line component 
is highly correlated with Wolf-Rayet (WR) wind-feature emission 
arising from a population of WNL stars, strongly suggesting the 
enhanced nitrogen to be the result of chemical enrichment by the 
N-rich winds of WN-type WR stars. This supports the postulate of 
nitrogen self-enrichment due to WR stars that had been suggested 
by Walsh & Roy (1987) and Pagel et al. (1986) to explain elevated 
N/O levels. 

James et al. (2010) also presented an integral field spectro- 
scopic study of the interacting dwarf galaxy UM 420 for which 
the long-slit spectroscopic study of Izotov & Thuan (1999) had re- 
ported log(N/0) = — 1 .08. By performing a spatially resolved anal- 
ysis of its physical and chemical conditions, James et al. (2010) 
found an average ratio applicable to the whole galaxy of log(N/0) 
= —1.64, i.e. 3.6 times lower than previous estimates. The VLT 
observations of James et al. (2010) did not confirm previous detec- 
tions of WR spectral features in UM 420, highlighting the difficul- 
ties involved in identifying WR stellar populations in distant dwarf 
galaxies. 

In this paper, we present IFU observations of another mem- 
ber of this subset of potentially N/O-anomalous galaxies, UM 448 
(UGC 06665, Mrk 1304). With a reported factor of ~3 nitrogen 
excess (Bergvall & Ostlin 2002), this galaxy presents an excel- 
lent opportunity to further our understanding of peculiar N/O ra- 
tios in some BCGs and the mechanisms responsible. UM 448 is 
a complex system, with a disturbed morphology possibly due to a 
merger event (Dopita et al. 2002). Thus the need for spatially re- 
solved spectroscopy is two-fold: firstly in determining the range of 
physical properties and chemical compositions that are present, and 
secondly, in spatially relating any chemical anomalies to environ- 
mental factors (e.g. WR populations). 

UM 448 is a starburst galaxy (Keel & van Soest 1992; Kew- 
ley et al. 2001) consisting of various star forming knots, diffuse 
Ha emission in the north, and a SE tidal extension (Dopita et al. 
(2002)). Broad He II A4686 emission was detected by Masegosa 
et al. (1991) and subsequently by Guseva et al. (2000), who re- 
ported an atypical weakness in this line in comparison with other 
WR-signature lines. As with Mrk 996 and NGC 5253, a WR popu- 
lation may well be the source of N-enhancement within this BCG. 

In this paper we present integral field spectroscopy (IFS) ob- 
servations obtained with the VLT-FLAMES IFU. Crucially, this in- 
strument is sufficiently blue-sensitive to detect the [O II] A3727 
doublet line, meaning that we are not required to adopt uncer- 
tain ionisation correction factors to determine oxygen abundances. 
These data afford us a new spatiokinematic 3-D view of UM 448, 
with the spatial and spectral resolution to undertake a full multi- 
velocity component analysis. We adopt a distance of 75.9 Mpc for 



UM 448 from the velocity measured here of +5579±5 km s 
(z = 0.018660 ± 1.6 x 10~ 5 ), with a Hubble constant of 
H =73.5 km s" 1 Mpc" 1 (DeBernardis et al. 2008). 



2 OBSERVATIONS AND DATA REDUCTION 

Observations of UM 448 were obtained with the Fibre Large Ar- 
ray Multi Element Spectrograph, FLAMES (Pasquini et al. 2002) 
at Kueyen, Telescope Unit 2 of the 8.2 m VLT at ESO's observatory 
on Paranal, in service-mode on dates specified in Table 1. Obser- 
vations were made with the Argus Integral Field Unit (IFU), with 
a field-of-view (FoV) of ll'.'5x7'.'3 and a sampling of O'.'52/lens. 
In addition to science fibres, Argus has 15 sky-dedicated fibres that 
simultaneously observe the sky and that were arranged to surround 
the IFU FoV at distances of ~35-58". The positioning of the IFU 
aperture for UM 448 is shown overlaid on SuSI2 images 1 in Fig. 1. 

Four different low-resolution (LR) gratings were utilised: 
LR1 (AA3620-4081, 24.9±0.3 kms" 1 FWHM resolution), 
LR2 (AA3964-4567, 24.7±0.1 kms" 1 ), LR3 (AA4501-5078, 
24.6±0.3 kms" 1 ) and LR6 (AA6438-7184, 21.9±0.4 kms" 1 ). 
This enabled us to cover all the important emission lines needed for 
an optical abundance analysis, as detailed in Section 3. The photo- 
metric conditions, airmass and exposure time for each data set are 
detailed in Table 1. In addition to the science frames, spectropho- 
tometric standard star observations (Table 1), as well as continuum 
and ThAr arc lamp exposures were obtained. 

The data were reduced using the girBLDRS pipeline (Blecha 
& Simond 2004) 2 in a process which included bias removal, lo- 
calisation of fibres on the flat-field exposures, extraction of indi- 
vidual fibres, wavelength calibration and rebinning of the ThAr 
lamp exposures, and the full processing of science frames which 
resulted in flat-fielded, wavelength-rebinned spectra. The pipeline 
recipes 'biasMast', 'locMast', 'wcalMast', and 'extract' were used 
(cf. Tsamis et al. 2008). The frames were then averaged using the 
'imcombine' task of IRAF 3 which also performed the cosmic ray 
rejection. The flux calibration was performed within IRAF using 
the tasks CALIBRATE and STANDARD. Spectra of the spectropho- 
tometric stars quoted in Table 1 were individually extracted with 
girBLDRS and the spaxels containing the stellar emission were 
summed up to form a ID spectrum. The sensitivity function was 
determined using SENSFUNC and this was subsequently applied 
to the combined UM 448 science exposures. The sky subtraction 
was performed by averaging the spectra recorded by the sky fi- 
bres and subtracting this spectrum from that of each spaxel in the 
IFU. Custom-made scripts were then used to convert the row by 
row stacked, processed CCD spectra to data cubes. This resulted in 
one science cube per grating, i.e. four science cubes for the source 
galaxy. 

1 SuSI2 images were obtained from the ESO archive, with total exposure 
times of 1200, 1200, 1800 and 2700 s, average airmasses of 1.34, 1.28, 1.22 
and 1.17 and image qualities of 0.75", 0.93", 0.87" and 0.56" FWHM for 
the R-, V-, B- and l/-band images, respectively. The average image qual- 
ity of the images is ~0.78" FWHM. WCS information was assigned via 
position matching with GSC2. Other than aligning and median-combining 
the separate frames using IRAF's IMALIGN and IMCOMBINE, respectively, 
no additional reduction was performed on the images. 

2 The girBLDRS pipeline is provided by the Geneva Observatory 
http://girbldrs.sourceforge.net 

3 IRAF is distributed by the National Optical Astronomy Observatory, 
which is operated by the Association of Universities for Research in As- 
tronomy 
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Figure 1. SuSI2 (Superb-Seeing Imager, installed on the ESO 3.5m New Technology Telescope (NTT) at La Silla Observatory, Chile) colour-composite 
images of UM 448 (red, i?-band; green, V-band; blue B-band), with the 11,5" X7.3" FLAMES IFU aperture overlaid. Each image has been scaled to 
enhance different features within UM 448. The left-panel shows two large knots of continuum or ionised gas emission that extend away from the galaxy in a 
south-westerly direction, outside the IFU FoV. The right-panel shows the two distinct portions of the galaxy, a blue eastern tail and a reddish knot of stellar 
continuum. SuSI2 images were obtained from the ESO archive, observed as part of programme ()78.B-0358(A). 



Table 1. FLAMES IFU observing log 

Date Grism Wavelength Range Exp. time Avg. Airmass FWHM seeing Standard Star 







(A) 


(s) 




(arcsec) 




24/05/2009 


LR1 


3620-4081 


3x233 


1.11 


0.85 


LTT7379 




LR2 


3964-4567 


3x300 


1.12 


0.87 


LTT7379 




LR3 


4501-5078 


3x230 


2.24 


0.80 


LTT7379 


18/06/2009 


LR1 


3620-4081 


3x233 


1.12 


0.47 


Feigel 10 




LR2 


3964-4567 


3x300 


1.13 


0.42 


Feigel 10 




LR3 


4501-5078 


3x230 


1.16 


0.50 


Feigel 10 


28/07/2009 


LR6 


6438-7184 


3x400 


1.75 


0.47 


EG21 


26/01/2010 


LR1 


3620-4081 


3x233 


1.15 


0.54 


LTT4364 




LR2 


3964-4567 


3x300 


1.13 


0.53 


LTT4364 




LR3 


4501-5078 


3x230 


1.11 


0.63 


LTT4364 


16/03/2010 


LR1 


3620-4081 


3x233 


1.10 


0.57 


LTT4364 




LR2 


3964-4567 


3x300 


1.11 


0.50 


LTT4364 




LR3 


4501-5078 


3x230 


1.12 


0.56 


LTT4364 


08/07/2010 


LR6 


6438-7184 


8x600 


1.63 


0.77 


Feigel 10 


09/07/2010 


LR1 


3620-4081 


3x233 


1.50 


0.77 


Feigel 10 




LR2 


3964-4567 


3x300 


1.60 


0.72 


Feigel 10 




LR3 


4501-5078 


3x230 


1.74 


0.66 


Feigel 10 



Observations of an object's spectrum through the Earth's at- 
mosphere are subject to refraction as a function of wavelength, 
known as differential atmospheric refraction (DAR). The direction 
of DAR is along the parallactic angle at which the observation is 
made. Following the method described in James et al. (2009), each 
reduced data cube was corrected for this effect using the algorithm 
created by Walsh & Roy (1990) ; this proceedure calculates the 
fractional pixel shifts for each monochromatic slice of the cube rel- 
ative to a fiducial wavelength (i.e. a strong emission line), shifts 
each slice with respect to the orientation of the slit on the sky and 
the parallactic angle and re-combines the DAR-corrected cube. All 
pixel shifts were < 1 pixel in size, indicating uncertainties due to 
extinction or temperature gradients are negligible. 



3 MAPPING OF LINE FLUXES AND KINEMATICS 



The distribution of H I Balmer line emission across UM 448 which 
is indicative of current massive star formation activity, and the mor- 
phology of the system as it appears on the SuSI2 images, were used 
to define three regions of roughly similar Balmer line luminosity. 
Our analysis of the galaxy is based on properties averaged over 
each of these three areas whose boundaries are displayed on the 
Ha map in Figure 2. Regions 1 and 2 are both contained in the blue 
eastern tail of UM 448 as seen on Fig. 1 (right), containing the bulk 
of Balmer line emission. Region 3 encompasses the interface of the 
blue eastern tail with the redder western part of the galaxy, and fur- 
ther contains the peak of the continuum emission in the vicinity of 
Ha (see Fig. 12). Whilst each of these regions shows a relatively 
small variation in surface brightness, the composite structure of the 
emission lines can vary considerably as illustrated by the Ha emis- 
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Figure 2. Ha emission per 0'.'52x(X'52 arcsec 2 spaxel. Three regions were defined across UM 448 over which summed spectra were extracted for analysis. 
On the right the Ho line profile from each summed region is shown, together with a two-component Gaussian fit. Blue dashed lines demarcate the wavelength 
range used for summing the flux. North is up and east is to the left. 





Fa 




[On] A3727 


98.02 ± 1.18 


125.40 ± 6.20 


[Oil] A3729 


135.65 ± 1.43 


173.55 ± 8.52 


[Nelll] A3868 


23.41 ± 1.71 


29.24 ± 2.54 


H8+Hel A3888 


15.47 ± 2.24 


19.26 ± 2.93 


HS 


21.57 ± 1.21 


25.70 ± 1.85 


H7 


43.33 ± 0.83 


48.96 ± 2.31 


[OIII] A4363 


1.90 ± 0.19 


2.14 ± 0.23 


He I A4471 


4.50 ± 0.49 


4.94 ± 0.57 


H/3 


100.00 ± 0.99 


100.00 ± 3.94 


[Olll] A4959 


91.86 ± 1.66 


89.76 ± 3.71 


[Nil] A6548 


15.83 ± 0.52 


11.66 ±0.49 


Ha 


407.05 ± 13.27 


301.47 ± 12.62 


[N II] A6584 


58.55 ± 0.51 


42.91 ± 1.17 


He 1 A6678 


2.86 ± 0.06 


2.07 ± 0.07 


[Sll] A6716 


49.17 ± 0.44 


35.38 ± 0.94 


[S ll] A6731 


36.07 ± 0.34 


25.95 ± 0.70 


c(H/3) 


0.42 ± 0.02 


F(H/3) 


24.48 ±0.17 



Table 2. UM 448 fluxes and de-reddenened line intensities (both relative 
to F(H/3)=I(H/3)=100) for summed spectra over all of UM 448. Line fluxes 
were extinction-corrected using the c(H/3) values shown at the bottom of 
the table, calculated from the relative Ha, H/3 and H7 fluxes. The line in- 
tensities listed here were used for summed-spectra T c and iV c diagnostics 
and regional ionic abundance calculations listed in Table 3. 



nents in the majority of emission lines seen in the 300 spectra of 
UM 448 across the IFU aperture. Following the methods outlined in 
James et al. (2009) and Westmoquette et al. (2007), we utilised an 
automated fitting proceedure called PAN (Peak ANalysis; Dimeo 
2005), that allowed us to read in multiple spectra simultaneously, 
interactively specify initial parameters for a spectral line fit and 
then sequentially process each spectrum, fitting Gaussian profiles 
accordingly. The output consists of the \ 2 value for the fit, the fit 
results for the continuum and each line's flux, centroid and FWHM, 
and the respective errors of these quantities which were propagated 
in all physical properties derived afterwards 4 . 

Even though some low S/N lines may in reality be composed 
of multiple emission components, it is often not statistically signif- 
icant to fit anything more than a single Gaussian. In order to rig- 
orously determine the optimum number of Gaussians required to 
fit each observed profile, we followed the likelihood ratio method 
outlined by Westmoquette et al. (201 1, their Appendix A). S/N ra- 
tio maps were first made for each emission line, taking the ratio 
of the integrated intensity of the line to that of the error-array pro- 
duced by the FLAMES pipeline, on a spaxel-by-spaxel basis. Sin- 
gle or multiple Gaussian profiles were then fitted to the emission 
lines, restricting the minimum FWHM to be the instrumental width 



sion line profiles displayed alongside the monochromatic line map 
shown in Figure 2. 

The high spectral resolution and signal-to-noise (S/N) ratio 
of the data allow the identification of multiple velocity compo- 



4 The uncertainties output by the pipeline were used as the starting point of 
the enor analysis; they also include uncertainties from the flat-fielding and 
the geometric correction of the spectra and not just photon-noise (Blecha & 
Simond 2004). 



(see Section 2). Suitable wavelength limits were defined for each 
line and continuum level fit. Further constraints were applied when 
fitting the [S II] and [O II] doublets: the wavelength difference 
between each line in the doublet was taken to be equal to the red- 
shifted laboratory value when fitting the velocity component, and 
their FWHMs were set equal to one another. 

The majority of high S/N ratio line profiles in the spectra 
of UM 448 are optimally fitted with a single narrow Gaussian 
(FWHM<100 kms , component CI hereafter), a single under- 
lying broad Gaussian component (FWHM 150-300 kms - , C2), 
and a third narrow Gaussian which can appear either red- or blue- 
shifted with respect to CI (C3 hereafter). In consideration of the 
different 'features' of each component (i.e. FWHM and/or veloc- 
ity), we consider each component as being real, i.e. arising from 
gas with different physical conditions. Towards the outer regions 
of the galaxy, where the S/N ratio decreases, the lines are optimally 
fitted with only CI and C2 velocity components. Low S/N ratio 
emission lines (e.g. H8+He I A3888) were optimally fitted with a 
single narrow Gaussian profile (CI). Figure 3 shows the complex 
profile of the Ha line, along with the spatial distribution of its con- 
stituent velocity components. 

Tables A1-A3 list measured FWHMs and observed and de- 
reddened fluxes for the fitted Gaussian components of the detected 
emission lines across the three regions of massive star formation 
defined in Fig. 2. The fluxes are for spectra summed over each re- 
gion and are quoted relative to the flux of the corresponding H/3 
component. They were corrected for reddening using the galactic 
reddening law of Howarth (1983) using c(H/3) values derived from 
the Ha/H/3 and H7/H/3 line ratios of their corresponding compo- 
nents, weighted in a 3:1 ratio °, respectively, in conjunction with 
theoretical Case B ratios from Storey & Hummer (1995). An ex- 
tinction map was also derived using the the Ha/H/3 and H7/H/3 
emission line ratio maps and is shown in Fig. 4. This was com- 
puted using the integrated line profiles (i.e. C1+C2+C3) to obtain 
a global extinction view of the galaxy. 

A Milky Way reddening of E(B - V) = 0.025 mag in the di- 
rection to UM 448 is indicated by the extinction maps of Schlegel 
et al. (1998), corresponding to c(H^)=0.036. Total c(H^) values ap- 
plicable to the galaxy and its individual emission line components 
in each main aperture-region are listed in their respective tables 
(Tables A1-A3). 

3.1 Ha maps 

The Ha surface brightness of UM 448 (Fig. 3a) displays a hor- 
izontal 'S'-shaped morphology, extending a total of ~ 10" east 
to west, almost symmetrically about the FoV center. The lower 
arm of the galaxy extends ~ 5" to the east and contains the 
galaxy's peak Ha emission (corresponding to designated region 
1). Based on our fitting method, the average surface brightness 
within region 1 can be de-composed to ~ 1.8 x 10" 14 and ~ 
2.8 x 10 -14 ergs cm -2 s -1 arcsec -2 for line components CI and 
C2 respectively. The upper arm is seen to extend ~ 5" in the oppo- 
site direction and displays an overall slightly lower surface bright- 
ness of ~ 1.2 x 10" 14 ergs cm" 2 s -1 arcsec -2 for both compo- 
nents. Component C3 appears in three separate regions; two regions 

5 An exception was made for region 3 where the low S/N ratio of H7 pro- 
hibited the robust decomposition into the three components observed in Ha 
and H/3. Therefore, only the Ha/H/3 line ratio was used for the calculation 
of c(H/3) for this region. 



A 3D Chemodynamical Study of UM 448 5 




-4 -2 2 4 
arcsecs 



0.00 0.19 0.38 0.56 0.75 0.94 
E(B-V) 

Figure 4. Map of E(B-V) for UM 448, derived from the integrated flux of 
H I emission lines. Integrated light Ho emission line contours are overlaid. 
North is up and east is to the left. 



either side of the central boundary region, i.e. between the two stel- 
lar components (see Fig. 1), and a third region along the most west- 
ern part of the upper arm. C3 had an average surface brightness of 
~ 4.4 x 10 -14 ergs cm" 2 s" 1 arcsec" 2 . The total integrated Ha 
flux for UM 448 is 9.24x 10 -13 ergs cm -2 s -1 , which is in good 
agreement with the value of ~ 11.5 x 10 -13 ergs cm -2 s -1 de- 
rived by Dopita et al. (2002) from narrow band Ha imaging. Rela- 
tive line intensities are listed in Tables A1-A3 for summed spectra 
over each of the three designated regions, for each separate velocity 
component. 

The Ha velocity map (Fig. 3b) displays a rather similar ra- 
dial velocity structure for both CI and C2 components. Both show 
a smooth velocity distribution between +50 and —180 kms -1 , 
which suggests solid body rotation within the galaxy. These maps 
can be used to derive an axis of rotation for the gaseous component 
of UM448 along position angle, PA ~ 36° (shown as a dashed line 
in Fig. 3b). Interestingly, the alignment of this axis also divides the 
blue eastern 'tail' from the red, western part of the galaxy (Fig. 1). 
The radial velocity structure of component C3 is less straightfor- 
ward than its counterparts. In the two regions near the center of 
UM 448, the velocities range from -100 to +100 kms -1 , whilst 
in the extreme west of the galaxy they are consistently negative 
down to —300 kms -1 . The velocity structure there could be the 
signature of a discrete blue-shifted outflow, which is perhaps re- 
lated to the detached knots of emission present in the SuSI2 images 
~ 10" from the southwest corner of the Argus aperture (Fig. 1) 
- i.e. we are detecting fast-moving gas, kinematically decoupled 
from the main body of UM448. 

The maps of the Ha FWHM show a slight spatial correlation 
with the flux distribution (Fig. 3c). The eastern regions displays 
decreased FWHMs for both the CI and C2 line components, in 
comparison to the upper arm structure: the average FWHM for the 
lower arm is -40-60 kms" 1 for CI and 120-160 kms" 1 for C2, 
which increases to ~90-100 km s" 1 and 200-240 km s" 1 , respec- 
tively, in the upper arm. Interestingly, the region of low FWHM 
is co-spatial with the main peak in flux within the eastern arm. If 
the overall appearance of UM448 is indicative of a close encounter 
or merger of two systems, then it would make sense for the gas 
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Figure 3. Maps of UM 448 in the Ho velocity components C1-C3: the top row shows an example 3-component fit to the Ha line and the structure of the 
velocity component shown in red is displayed below; a) flux in units of 10 — 16 ergs s~ 1 cm -2 arcsec~ 2 per 0'.'52x0'. / 52 spaxel; b) radial velocity relative 
to the heliocentric systemic velocity of +5579 kms —1 (the dashed line refers to the PA~36° rotation axis discussed in the text; c) FWHM corrected for the 
instrumental PSF. See text for details. Maps are shown for S/N > 10 in the integrated line profile. North is up and east is to the left. 



turbulence, and hence the FWHM, to increase as one approaches 
the transition region between the two colliding components (the 
blue/red interface in Fig. 1 right). On the other hand, component 
C3 shows a typical FWHM of 70 km s _1 suggesting an association 
with more orderly gas volumes. 

3.2 Forbidden-line maps 

Emission line maps of UM 448 in the light of the [O III] [N II] 
and [Ne III] lines for the various velocity components are shown 
in Fig. 5. The morphology of [O III] A4959 is similar to that of 
Ha, with the surface brightness of the lower arm being higher by a 
factor of ~2 than that of the upper arm (for component CI). In con- 
trast, [N II] displays a rather constant surface brightness throughout 
either arm, with a peak in emission within our designated region 2. 
For both [O III] and [N II] line component C2 displays a simi- 
lar morphology to the CI counterparts, as was the case with Ha. 



Component C3 of [O III] A4959 has a wider spatial extent than 
its counterpart in Ha. For [N II] component C3 is seen to overlap 
roughly with the central and eastern regions observed for Ha-C3. 
In the central regions C3 shows symmetry approximately about the 
rotation axis defined in Fig 3b; this may a signature of an outflow, 
or it simply reflects the highly disturbed kinematics in the "collision 
zone" of the merging galactic components at the Region l/Region2 
aperture interface. [O III] A4363 has also been detected in veloc- 
ity components CI and C2 whose peak surface brightness is found 
within region 3 (see Fig. 10 for an example fit to the summed spec- 
tra across this region), whereas A4959 peaks in region 1. This is 
mirrored in the high electron temperature of region 3 versus region 
1 (see below). The [Ne III] A3868 emission also peaks within re- 
gion 1, with more diffuse emission extending throughout regions 2 
and 3. 

UM 448 is shown in the light of the [O II] and [S II] dou- 
blet lines in Fig. 6. Due to their relatively low S/N ratio (compared 
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Figure 5. Dereddened intensity maps of UM 448 in [O III] A4959, [N II] A6584, [O III] A4363 and [Ne III] A3868 for the separate velocity components CI, 
C2 and C3. [O III] A4363 is shown after being re-binned by 1,5x1.5 spaxels and re-mapped onto the original grid (see the text for details). Only spaxels with 
S/N> 10 are shown for [O III] A4959 and [N II] A6584, >5 for [Ne III] A3868 and > 1 for [O III] A4363. North is up and east is to the left. 



to the H I Balmer lines and [O III] A4959), both doublets could be 
decomposed to only the CI and C2 velocity components (Fig. 7). 
The [O II] doublet, with a rest-wavelength separation of 2.783 A 
was particularly hard to decompose into a narrow+broad compo- 
nent since the large FWHM of the component C2 would effectively 
merge the doublet into a single line. To alleviate this, each spaxel 
with a broad component [O II] line detection was examined and 
only selected when there was sufficient S/N ratio to allow for a de- 
blended fit (see Fig. 7a). The larger wavelength separation of the 
[S II] doublet allowed the broad component (C2) to be also fit- 
ted in relatively low S/N ratio spaxels (e.g. Fig. 7b), and was thus 
mapped more extensively than [O II] (compare the lower panels 
of Fig. 6a, b). The [O II] A3727+A3729 emission has a rather dis- 
jointed morphology compared to the integrated Ha surface bright- 
ness. The ratio of the [O II] doublet components commonly used 
as a density diagnostic, (Fig. 6, left panels) displays negligible 
structure in relation to the overlaid Ha. In contrast to [O II], the 
[S II] A6716+A6731 doublet has a morphology that is more akin 
to Ha with three distinct peaks in flux. The A6716/A6731 ratio 
(Fig. 6, right panels), which is also used as a density diagnostic, 
decreases in regions that align with the peaks in Ha emission, and 
this is indicative of increased electron density there. 



4 ELECTRON TEMPERATURE AND DENSITY 

The [O III] (A5007 + A4959)/A4363 intensity ratios were used to 
determine electron temperatures. Two emission line doublets are 
available to determine the electron density, N a ; namely the [O II] 
A3726/A3729 and [S II] A6716/A6731 ratios. T c and N c values 
were computed using IRAF'S ZONES task in the NEBULA pack- 
age. This makes use of technique that derives both the tempera- 
ture and the density by making simultaneous use of temperature- 
and density-sensitive line ratios. Atomic transition probabilities for 
2+ andO+ were taken from Wiese et al. (1996) whilst collision 
strengths were taken from Lennon & Burke (1994) and McLaugh- 
lin & Bell (1993), respectively . For [S II] we used the transitional 
probabilities of Keenan et al. (1993) and the collisional strengths of 
Ramsbottom et al. (1996). 

Where possible T e and N e values were derived for each ve- 
locity component within each of the aperture regions defined in 
Fig. 2 and are listed in Table 3. The mean electron densities from 
the [S II] and [O II] ratios agree to within ~20%; the [O II] density 
map is however significantly more noisy. The [S II] densities were 
therefore adopted for the analysis. 

T e and N a maps computed for line components CI and C2 are 
shown in Fig. 8, and mean values for each of the aperture regions 
1-3 of Fig. 2 are listed in Table 3. Fig. 8 (left) shows evidence for 
a small amount of T c variation throughout UM 448. Aperture re- 
gions 1 and 2 show temperatures of ~ 12,000-13,000 K. Region 3 
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Figure 6. Dereddened emission line maps of UM 448 in doublet ratioes [O II] A3727/A3729 (top-panel) and [O II] A6716/A6731(bottom-panel), in the separate 
Cl(narrow) and C2 (broad) velocity components, with integrated Ha emission line flux overlaid in contours. North is up and east is to the left. 



shows temperatures of 15,300 K for the narrow line component 
CI and 10,200 K for the broad line component C2. iV B variations 
are also apparent and denser areas are correlated with peaks in Ha 
emission. Regions 1 and 2 have an average 7Y C ~120-150 cm -3 , 
whereas region 3 is only marginally denser with 200 cm -3 . 



5 CHEMICAL ABUNDANCES 

Ionic abundance maps relative to H + were created for the N + , + , 
2+ , Ne 2+ and S + ions, using the AA6584, 3727+3729, 4959, 
3868 and 67 17+673 1 lines respectively. Examples of the flux maps 
used in the derivation of these ionic abundance maps can be seen 
in Figure 5. Abundances were calculated using the ZONES task 
in IRAF, using the respective T and 7V C maps described above, 
with each FLAMES spaxel treated as a distinct 'nebular zone' with 
its own set of physical conditions. These results were compared 
with those obtained using the multilevel ionic equilibrium program 
EQUIB (originally written by I. D. Howarth and S. Adams at Univ. 
College London), which comprises a different set of atomic data 
tables. Typical differences in the abundance ratios were found to be 
5% or less. To check the validity of using T c ([0 III]) for the singu- 
larly ionized elements (N + , S + , and + ), we also computed abun- 
dance maps useing a T c ([0 II]) map (created using the relationship 
with T c ([0 III]) described in Izotov et al. (2006)). Discrepancies 
between the two temperatures were found to be low (~ 100-400 K) 
and typical differences in ionic abundance were ~2-7%, These 



errors are lower than those caused by uncertainties in T c ([0 III]) 
alone and we therefore do not include it as a source of error. 

Ionic nitrogen, neon and sulphur abundances were converted 
into N/H, Ne/H and S/H abundances using ionisation correction 
factors (ICFs) from Kingsburgh & Barlow (1994).The O/H abun- 
dance was obtained by adding the 2+ /H + and + /H + abundance 
maps (i.e. assuming that the O /H + ratio is negligible). Since the 
[S III] A6312 line was not within the wavelength range of our data, 
the S 2+ /H + abundance was estimated using the empirical relation- 
ship between the S 2+ and S + ionic fractions from the corrected 
equation A38 of Kingsburgh & Barlow (1994), namely S 2+ /S + 
= 4.677 x (0 2+ /0+) 433 . The S/H results should therefore be 
considered with caution, but they are not at variance with those 
obtained using the ICF prescriptions of Izotov et al. (2006). O/H, 
N/H and log(N/0) abundance maps are shown in Fig. 9 and are 
described below. Regional abundances derived directly from tak- 
ing error-weighted averages over the abundance maps are listed in 
Table 3. 

5.1 Maps 

The elemental oxygen abundance map shown in the middle panel of 
Figure 9 displays a single peak whose location correlates with the 
peak in Ha relating to region 1 . Minimal abundance variations are 
seen across UM 448, with a slightly decreasing oxygen abundance 
as one moves from regions 1-3, with an overall average abundance 
for CI of 12 + log(0/H) = 8.15 ± 0.10. Adopting a solar oxy- 
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Figure 8. [O III] electron temperature maps of UM 448 for CI (narrow) and C2 (broad) velocity components. Overlaid circles represent the size of the T a 
uncertainty within each spaxel (percentage errors are scaled by the 0.52" diameter of a spaxel). North is up and east is to the left. 
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Figure 7. Example single-spaxel spectra and fits to the 
[O II] AA3727,3729 and [S II] AA67 16,6731 lines. Flux is in units of 
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ergs cm 2 s 1 arcsec 2 A 1 and wavelengths are in the 



rest-frame. 



gen abundance of 8.7 1±0. 1 relative to hydrogen (Scott et al. 2009), 
this corresponds to ~0.28 Zq. This is comparable to the value ob- 
tained by Izotov & Thuan (1998, IT98 hereafter), who measured 
an abundance of 12 + log(0/H) = 7.99±0.04. The C2 veloc- 
ity component displays a higher oxygen abundance than the Cl- 



component gas within region 3 by ~0.4 dex, which, judging from 
the agreement between CI and C2 [O III] 4959 A flux (Fig. 5) is 
mostly attributable to the T c being ~5000 K cooler (if the CI re- 
gion was instead much denser but of similar oxygen abundance to 
the gas sampled by the C2 line component one might expect the 
A4959/H/3 CI ratio to be much smaller something which is not ob- 
served - Table A3). 

The upper panel of Fig. 9 shows the N/H abundance across 
UM 448. The CI peak in nitrogen abundance is offset from that of 
oxygen, peaking at the location of region 2 with an average N/H 
ratio of 1.04±0.01 x 10~ 5 . For region 1 it is slightly lower at 
9.40±0.02 x 10~ 6 and for region 3 even lower at 7.00±0.03 x 
10~ 6 . The C2 component gas in region 3 has a nitrogen abundance 
that is higher than its CI counterpart by a factor of 2. As a result, 
there is a peak in the N/O ratio within UM 448 located at region 2, 
which can be seen in the bottom panel of Fig. 9. For CI, regions 1 
and 3 agree with an average of log(N/0) ~ -1.23±0.04, whereas 
region 2 has an average (N/O) ratio of ~ -1.03±0.05. The map of 
log(N/0) shows that this peak actually reaches ~ —0.80 at the pre- 
cise location of region 2 (as can be seen from the Ha contours). 
The N/O ratio of the (broad) C2 component gas in region 3 (~ - 
1.26±0.13) is in good agreement with the average N/O of CI gas in 
that region. The long-slit study by IT98 measured an N/H ratio of 
9.5 ±2.0 x 10 -6 , in agreement with the averages of regions 1 and 2 
derived from the N/H map shown in Fig. 9. The average CI N/O ra- 
tios for region 3 is slightly higher (by ~0.2 dex) than those of other 
galaxies of similar oxygen metallicity (see Lopez-Sanchez & Es- 
teban 2010, their figure 11). They are, however, consistent within 
the uncertainties. Region 2, however, has a higher value than the 
average value of ~ —1.4 ±0.2 expected for its metallicity, even 
when taking the uncertainties into consideration, for both the av- 
eraged value of region 2 and on a spaxel-by-spaxel basis. For this 
region we do observe the factor of ~ 3 nitrogen excess reported by 
Pustilnik et al. (2004). 

The Ne/H and S/H abundance distribution throughout UM 448 
for the narrow (CI) emission line component show only minimal 
abundance variations in relation to the Ha emission. Ne/H values 
averaged over each region range from ~2.06-6.55 x 10 -5 , with the 
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Figure 9. UM 448 nitrogen, oxygen and log(N/0) abundance maps for velocity components CI and C2, with integrated light Ha emission line contours 
overlaid in dashed line. North is up and east is to the left. 



highest Ne/H abundance aligning with region 2. The Ne/O ratios 
of region 1 and 3 are typical for those expected by BCGs at simi- 
lar metallicities; the averaged region values lie within the reported 
range of log(Ne/0) ~ -0.9 to -0.7 (Izotov et al. 2006). However, 
the Ne/O ratio for region 2 is ~0.5 dex higher than the expected 
Ne/O value. Calculations of Ne/H were also made using the Ne 2+ 
ICF of Izotov et al. (2006), however the Ne/O ratio for region 2 
remained ~0.2 dex higher than the expected value. The S/H distri- 
bution across UM 448 is relatively constant for the CI and C2 gas 
components, with each region displaying averages within the range 
of ~ 4.7 to 7.4 x 10~ 6 . The CI S/O ratio for all regions are slightly 
higher than expected for blue compact galaxies, log(S/0)~ —1.9 



to —1.5 (Izotov et al. 2006), and the S/O value C2 component in 
Region 3 lies at the upper limit of this average. 



5.2 Comparison of spatially resolved and global spectra 
abundances 

Also listed in Table 3 are the T a , N c , ionic and elemental abun- 
dances derived from summed spectra over the entire galaxy (using 
line fluxes listed in Table 2). This particular set of results allows us 
to assess if global galaxy spectra can reliably represent the physi- 
cal properties of its ISM. For example, values derived from long- 
slit observations may contain luminosity-weighted fluxes and con- 
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Figure 10. [O III] A4363 profile observed in the spectra summed over re- 
gion 3, displaying a narrow (CI) plus broad (C2) component Gaussian pro- 
file. Flux, FWHM and velocity parameters for each component are listed in 
Table A3 



sequently erroneous ratios, or large apertures may include a mix- 
ture of gas with different ionization conditions and/or metallicities. 
This question was addressed by Kobulnicky et al. (1999) to assess 
whether global emission-line spectra can be depended upon to reli- 
ably recover the chemical properties of distant star-forming galax- 
ies. Overall, they find that when using direct-method techniques, 
global spectra can give an accurate representation of the oxygen 
abundances derived from individual H II regions, but only after ap- 
plying a correction of A(OfH)— 0.1 to the global spectra results. 
In addition to this, Pilyugin et al. (2012) found that oxygen abun- 
dances were systematically underestimated (and N/O consequently 
over-estimated) when using global SDSS (Sloan Digital Sky Sur- 
vey) spectra of H II regions with different physical properties. 

For UM 448, the temperature derived from the integrated spec- 
trum is in agreement with the flux-weighted mean T c across the 
three regions (14010 ± 2100 K), suggesting that the integrated 
[O III] ratio may have been more heavily weighted by areas of 
stronger [O III] A4363 emission. The oxygen abundance from the 
integrated spectrum, 12+log(O/H)=8.34±0.07, is higher than all 
CI -component values, suggesting that it may be weighted towards 
the high O/H abundance seen in the C2-component of region 2. Be- 
ing +0.18 dex higher than the flux-weighted average of the three 
regions (12+log(O/H)=8.16±0.18), the globally-derived (i.e. de- 
rived from the integrated spectrum) O/H is not fully representa- 
tive of the emission from the whole galaxy, even after applying the 
correction suggested by Kobulnicky et al. (1999). Nitrogen too, dis- 
plays a higher abundance when derived from the summed spectrum 
rather than the regional averages, and results in a (log) nitrogen-to- 
oxygen ratio of —1.15, which is approximately the expected value 
at this metallicity. The Ne/H abundance ratio from the integrated 
spectrum is almost a factor of 2 higher than the flux-weighted aver- 
age of the three regions (~2.14x 10 -5 ), whilst the log(Ne/0) ratio 
is in good agreement. The same applies for the S/H abundance, al- 
though less severe, where the globally-derived S/H ratio is a factor 
~1.5 higher than the flux-weighted average S/H across the three 
regions, but the log(S/0) value is in agreement between the two 
methods. 

Comparisons can also be made with the results from a long- 
slit study by IT98, which utilized summed spectra over a 2" x 7" 
region 'across the brightest region of the galaxy'. From this spec- 
trum they derive an T e of 12100±500 K (in line with the CI com- 
ponent values of our region 1), O, N, Ne and S, abundances that are 
systematically lower than those derived from our summed spec- 
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Figure 11. NTT-SuSI2 combined R, V and S-band images of UM 448 
with contours showing the continuum near Ha (top panel) and integreated 
Ho emission (bottom panel; see Fig. 3 for maps of the separated emis- 
sion line components). North is up and east is to the left. Ha contours are 
shown for the range 35.0 to 128.6 in steps of 22.0, in flux units of 10~ lfi 
ergs s — 1 cm -2 arcsec - 2 . The continuum peak region is represented by 
the white-dashed box; spectra were summed over this region to investigate 
stellar absorption features and averages in EW(H/3). 



trum and abundances relative to oxygen that are in relatively good 
agreement within the uncertainties. However, without knowing the 
exact positioning of this slit, a truly meaningful comparison cannot 
be made. 

Overall it can be said that whilst abundance ratios of heavy el- 
ements can be reasonably well obtained by summed spectra and/or 
long-slit observations, the ionic and elemental abundances relative 
to hydrogen are not well reproduced; this may be because the for- 
mer properties do not have a very severe sensitivity to T c biases, 
while the latter do. This is in addition to the fact that, in contrast 
to integral field spectroscopy, global properties derived from light- 
averaged methods cannot resolve abundance variations and any po- 
tential sites of enrichment are being averaged out. Of course, un- 
certainties will still exist within IFU observations due to averaging 
along the line of sight of individual spaxels; this is something that 
can only be remedied with the development of higher spectral res- 
olution IFUs. 



6 STELLAR POPULATIONS 
6.1 Age of the Stellar Population 

Figure 1 1 shows a combined NTT-SuSI2 R, V and B-band image 
that shows the stellar continuum within UM 448. Overlaid on the 
top panel are the contours of the continuum emission in the Ha 
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Table 3. Ionic and elemental abundances for UM 448, with the corresponding T e and 7V e used in their derivation. Columns 2-5 refer to abundances derived 
from maps, weight-averaged over regions defined in Fig. 2. Column 6 refers to T e , N e and abundances derived from spectra integrated over UM 448 
(corresponding fluxes are listed in Table 2). All ionisation correction factors (ICFs) are taken from Kingsburgh & Barlow (1994). S/H ratios were calculated 
using the ICF(S) values listed along with a predicted value of S 2 + (see text for details). 



spectral vicinity (~6 A windows either side of the line) that are 
in good agreement with the morphology of the broad-band images. 
In the lower panel we show the same broad-band images with the 
Ha integrated-light emission line contours. Significant differences 
between the morphology of the ionised gas emission and that of 
the stellar continuum are clearly apparent, especially in the top, 
western 'arm'. A strong continuum peak can be seen, which is not 
mirrored by the Ha emission, suggesting that a strong stellar com- 
ponent exists within this region. From the broad-band colour im- 
age shown in Figure 1 it can be seen that this peak in continuum is 
red in colour, signifying an older stellar population than that seen 
in the bluer, eastern 'arm'. Such a distinct, two component stellar 
structure, in conjunction with the difference in radial velocity be- 
tween them (see Fig. 3b), supports the interpretation that UM 448 
is a merging system comprised of two galaxies, with very differ- 
ent ages. Spectra summed over this continuum peak region (rep- 
resented by the white-dashed box overlaid in Fig. 1 1 -lower panel) 
do not reveal any stellar absorption features around the Ha, H/3, 
H 7 or HS lines, however Ca II H, K (AA3970, 3934) interstellar 
absorption features were detected. 

The age of a young population in a galaxy can be estimated us- 
ing hydrogen recombination lines, since they provide an estimate of 
the ionizing flux present, assuming a radiation-bounded HII region 
(Schaerer & Vacca 1998). In particular, the equivalent-width (EW) 
of H/3 can be used as an age indicator of the ionizing stellar popu- 
lation at a given metallicity. A map of EW(HP) for UM 448 can be 
seen in Figure 12. Whilst there are no clear peaks in EW(Hf3), the 
main region of decreased EW correlates with the largest peak in 
continuum emission (red contours) and that corresponds to region 
3. 

Following the method outlined by James et al. (2010), we can 



use this map, in conjunction with the metallicity map described 
in Section 5, to estimate the age of the most recent star forming 
episodes throughout UM 448, by comparing the regional observed 
average EW(H/3) values with those predicted by the spectral syn- 
thesis code STARBURST99 (Leitherer et al. 1999). For the mod- 
els we chose a metallicity 0.2 Z Q , the closest representative to the 
weighted average metallicity of UM 448 (~0.25Z Q , derived from 
regional averages), together with assumptions of an instantaneous 
burst with a Salpeter initial mass function (IMF), a total starburst 
mass of 1 X 10 6 M Q , a 100 M stellar upper mass limit (which ap- 
proximates the classical Salpeter (1955) IMF). Models were run for 
Geneva tracks with "high" mass loss rates and Padova tracks with 
thermally pulsing AGB stars included. For each of these evolution- 
ary tracks, two types of model atmosphere were used; firstly the 
Pauldrach-Hillier (PH) and secondly Lejeune-Schmutz (LS). The 
latter was chosen because it incorporates stars with strong winds, 
which would be representative of the WR population within the 
galaxy (if present). The stellar ages predicted by each model and 
the observed average EW(H/3) within each peak emission region 
in UM 448 were derived from Figure 13 and are listed in Table 4. 

The difference between the ages predicted by the Geneva 
and Padova stellar evolutionary tracks is relatively small, with the 
Geneva tracks predicting lower ages by up to 23%. The differ- 
ence in ages predicted by the PH and LS atmospheres are smaller 
still, with LS predicting lower ages by ~5%. We cannot com- 
ment on which atmosphere model is more appropriate given that 
the existence of WR stars within UM 448 is not confirmed by the 
current study, hence we adopt average stellar ages from the four 
model combinations. The current (instantaneous) star formation 
rates (SFRs) based on the Ha luminosities were calculated fol- 
lowing Kennicutt (1998) and are given in Table 4. SFRs corrected 
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EW (H/S) / A 

Figure 12. Map of the equivalent width of H/3 accross UM 448. Overlaid 
red lines are contours of H/3-region continuum; black lines are H/3 flux 
contours from Fig. 2. 



Figure 13. EW(Hf3) as a function of age, as predicted by the STAR- 
BURST ^ code for metallicities of 0.2Zq using a combination of Geneva or 
Padova stellar evolutionary tracks and Lejeune-Schmutz (LS) or Pauldrach- 
Hillier (PH) model atmospheres. The observed average J?W(H/))'s for each 
star forming region in UM 448, along with the peak in continuum high- 
lighted in Fig. 1 1 (lower panel), are overlaid (solid line). 



for the sub-solar metallicities are also given, derived following the 
methods outlined by Bicker & Fritze-v. Alvensleben (2005). 

Regions 1-3 contain ionizing stellar populations with aver- 
age ages of 4.6, 5.5 and 6.0 Myr, respectively, and the region of 
peaked continuum (Fig. 11) has an average age of 8.2 Myrs (Ta- 
ble 4). since Region 1 shows the highest ionizing flux, we would 
expect it to contain the youngest stellar population. It is notable 
that the ages of the stellar population increase through Regions 1-3 
as one moves away from the blue continuum peak (Fig. 1). How- 
ever, the SFR does not show a smooth decrease in the same direc- 
tion, decreasing slightly for region 2 before increasing again for 
region 3. Interestingly, and not surprisingly, the main red-colored 
stellar continuum peak (Figs. 1 and 11) displays an increase in age 
of almost 4 Myrs from the blue continuum peak (region 1) and a 
much decreased SFR compared to the three main regions of star- 
formation. Using narrow-band Ha imaging, Dopita et al. (2002) 
derived a total SFR of 6.6 M©/yr, which is in good agreement with 
that derived from our total Ha flux, 5.0 M0/yr. Accounting for its 
~0.2Zotnetallicity, this translates to 2.8 M©/yr. The evidence sug- 
gests that UM 448 is made up of two separate bodies undergoing a 
merger; the first is sampled by our aperture regions 1 and 2, where 
relatively recent bursts of star-formation have occurred. The second 
body with which it is merging may be an older galaxy (represented 
by the red continuum peak) sampled by aperture region 3, within 
which the collision could have triggered star-formation, hence an 
increased SFR but overall age that is representative of an aggregate 
old+young population. 

6.2 WR stars and N-enrichment 

Contrary to the conclusions of previous long-slit studies (Masegosa 
et al. 1991; Schaerer et al. 1999, and references therein) and in 
agreement with Guseva et al. (2000), our FLAMES IFU spectra 
do not reveal evidence of 'strong' WR emission bands at ~4640- 
4686 A. Unfortunately the C IV 5801,5812 A feature is not included 
in our wavelength coverage, although these features are usually 
significantly weaker than the ~4686 A feature. We find no evi- 
dence of a WR population in our single spaxel spectra and only 
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Figure 14. Sections of continuum-subtracted FLAMES IFU spectra, show- 
ing the blue 'WR bump' region within summed spectra over the entire 
galaxy and separate aperture regions 1, 2 and 3 (see Fig. 2). There are min- 
imal detections of WR emission features within both the spectra summed 
over the entire galaxy and also each of the individual aperture region spec- 
tra. Overlaid are the estimates of 4000±1100, 900±600, 300±600 and 
1300±400 WN (5-6) stars (solid red line), for the entire galaxy and regions 
1—3, respectively, using the templates of Crowther & Hadfield (2006), along 
with the Ict errors (dashed red line). Vertical dashed lines indicate the loca- 
tion of [Fe III] 4658 A and a weak, narrow, He II emission line at 4686 A. 
The dashed horizontal lines in the plots show the zero flux levels. 

minimal evidence in our spectra summed over the entire galaxy. 
When searching for WR features one must consider the width of the 
extraction aperture, which can sometimes be too large and dilute 
weak WR features by the continuum flux (Lopez-Sanchez & Este- 
ban 2008). This has indeed been the case for several other BCGs, 
where previous reports of detections of WR features in long-slit 
spectra were not confirmed by IFU spectra (e.g., James et al. 2010; 
Perez-Montero et al. 201 1). We therefore also examined the spectra 
summed over each individual region (Figure 2), which are shown 
in Fig. 14. 

In order to determine the number of WR stars from the 
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Table 4. Age of the latest star formation episode and current star formation rates. 







Starburst aj 


jes (Myr) 




Model 


Reg. 1 


Reg. 2 


Reg. 3 


Cont m Peak 


Padova-AGB PH 


5.06± 0.24 


5.86+ 0.13 


6.36+ 0.48 


8.61+ 0.74 


Padova-AGB LH 


4.96± 0.24 


5.66+ 0.13 


6.01+ 0.33 


8.51± 0.79 


Geneva-High PH 


4.41+ 0.27 


5.31± 0.16 


5.91+ 0.56 


7.91+ 0.74 


Geneva-High LS 


3.86+ 0.48 


5.11± 0.16 


5.61+ 0.50 


7.76+ 0.82 






Star formation rates (Mq yr 




SFR(Ho) a 


1.49±0.03 


0.75+0.02 


0.90±0.02 


0.34±0.01 


SFR(Ha) 6 


0.82±0.02 


0.41±0.01 


0.50+0.01 


0.18±0.01 



a Derived using the relationship between SFR and L(Ha) from Kennicutt (1998). 

b Corrected for sub-solar metallicities following Bicker & Fritze-v. Alvensleben (2005). 



strength of the 4640-4690 A feature for the summed spectra over 
regions 1-3, the local continuum was subtracted. A linear interpo- 
lation to the continuum from both sides of the 4640-4686 A region 
was employed based on the mean flux in a blue region (rest wave- 
lengths 4609 - 4640 A, avoiding the [Fe III] emission line and pos- 
sible C III absorption at 4647 A) and a red region (4711-4742 A). 
Utilizing an LMC WR (WN5-6) spectral template from Crowther 
& Hadfield (2006), the best match of the scaled WR star flux to 
the flux over the observed broad 4686 A region was determined by 
minimizing \ 2 over the peak of the WR feature, typically 4670 - 
4900 A. Since errors were available on the flux points these were 
used to generate a series of Monte Carlo realizations of each re- 
gion spectrum and the WR template flux was best fitted to each. 
The distribution the WR template flux to each realization was fitted 
by a Gaussian enabling the ler errors to be estimated. The flux of 
the WR template was then converted into the number of WR stars 
allowing for the distance of UM 448. We estimate 900 ± 600, 300 
± 600 and 1300 ± 400 WN stars for regions 1-3 respectively. A 
total population of 4000 ± 1 100 WN stars was estimated applying 
the same procedure for the summed spectrum for the entire galaxy. 
These error estimates are based purely on the flux errors and do not 
take into account any systematic error arising from the, necessarily 
simplified, removal of a linear continuum under the WR line.The 
WR template fits and the error ranges for the summed spectra over 
the entire galaxy and the three regions 1-3 are shown in Fig. 14. 
However, considering the large uncertainties, the existence of a WR 
population within UM 448 is still considered to be speculative and 
as a result, we cannot conclude that the Alog(N/0)~ +0.6 dex ex- 
cess observed in region 2 is due to nitrogen enrichment from WR 
stars. 

However, following the methodology outlined in Perez- 
Montero et al. (201 1), we can use the size of the N-enriched region 
in UM 448 (~1.6" x 1.6" = 0.6 kpc 2 ) in region 2 as a scale length 
to assess whether the reported WR stars could be the actual source 
of this N-pollution. Using modeled chemical yields at ~O.2Z , 
those authors study the effects of elements ejected into the ISM via 
winds as a function of time and pollution radii (see Perez-Montero 
et al. 201 1, section 4.3 for a full description of the model parame- 
ters). After ~4 Myr, the appearance of WR star winds increases the 
N/O ratio, where the value of N/O ratio only reaches that of region 
2 at radii much lower (e.g., r = 10 or 100 pc) than those mea- 
sured from our IFU observations. Therefore, although WR stars 
have been detected in previous long-slit studies, and are tentatively 
detected here, the high N/O ratio observed within region 2 is not 
likely produced by the pollution from stellar wind ejecta coming 



from WR stars. Instead, this enrichment may be related to other 
global processes occurring within UM 448, namely the likely in- 
teraction and/or merger between two bodies as inferred from the 
SuSI2 images and radial velocity maps. Studies suggest that in- 
teracting galaxies fall >0.2 dex below the mass-metallicity rela- 
tion of normal galaxies due to tidally induced large-scale inflow of 
metal-poor gas to the galaxies' central regions (e.g. Kewley et al. 
2006; Michel-Dansac et al. 2008; Peeples et al. 2009; Rupke et al. 
2010; Torrey et al. 2012). In addition, models by Koppen & Hensler 
(2005) have shown that a rapid decrease in oxygen abundance dur- 
ing an episode of massive and rapid accretion of metal-poor gas 
is followed by a slower evolution which leads to the closed-box 
relation, thus forming a loop in the N/O-O/H diagram. Such a sce- 
nario, in conjunction with metal-rich gas loss driven by SN winds, 
has been suggested by Amorin et al. (2010) to explain the system- 
atically large N/O ratio found within green pea galaxies. 



7 SUMMARY & CONCLUSIONS 

We have analysed FLAMES-IFU integral field spectroscopy of the 
BCG UM 448 and studied its morphology by creating monochro- 
matic emission line maps. This system shows signs of interaction 
and/or perturbation and it is possibly undergoing a merger. Based 
on the emission line morphology, we defined three regions of high 
Ha surface brightness which we use to derive error-weighted av- 
erages of the physical and chemical conditions across UM 448. 
Region 1 is to the east and corresponds to the tidal tail, region 3 
extends towards to the upper-west and region 2 lies between them. 
UM 448 exhibits complex emission line profiles, with most lines 
consisting of a narrow, central component and an underlying broad 
component, and a third, narrow blue-shifted component is also de- 
tected in Regions 2 and 3. 

Abundance ratio maps of O/H, N/H, Ne/H and S/H were cre- 
ated based on the direct method of estimating electron tempera- 
tures from the [O III] A4363/A4959 line ratio. Minimal variance in 
oxygen abundance are seen across UM 448 for the narrow velocity 
component, with an average abundance of 12+log(O/H)=8.15±0.1. 
However, in an area of UM 448 where the [O III] 4363 A line is re- 
solved into a narrow+broad component, an ~0.4 dex increase in 
oxygen abundance attributable to a decrease in T c of ~5000 K is 
found in the broad line region. We detect the previously reported 
high N/O ratio, with Alog(N/0)~ +0.4, in our aperture region 
2 only, spatially coincident with a peak in nitrogen abundance. 
Abundances derived from an integrated spectrum of UM 448 are 
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not a reliable representation of the ionic and elemental abundances 
throughout UM 448, which may be a consequence of the variance 
in T c throughout the nebular emission line components of UM 448 
detected in these integral field spectra. However, abundances rela- 
tive to oxygen from the two methods (i.e. map-based averages vs. 
integrated spectra) are consistent within the uncertainties. These 
results are significant when assessing the reliability of abundances 
derived from long-slit and integrated spectra of both nearby and 
high-redshift star-forming galaxies that do not contain homoge- 
neous ionization conditions and/or metallicities. 

Continuum images of UM 448 display a rather different mor- 
phology to those of the emission lines, with a peak in stellar light 
situated in between the second and third custom-aperture regions. 
The colour of this continuum peak (from SuSI2 broad-band im- 
ages) shows signs of a red, older stellar population in contrast to 
a younger, bluer stellar population situated in the eastern body of 
UM 448. Star formation rates were computed from the Ha line lu- 
minosities within each sampled star forming region and the ages of 
the last major star formation episode were estimated by fitting the 
observed Balmer line equivalent widths with STARBURST99 mod- 
els. The two merging sections of UM 448 have SFRs that differ by a 
factor of ~5 and starburst episodes separated by 4 Myr, providing 
evidence of a distinct age gradient between the two components. 
The latest major star formation event took place ~4.5 Myr ago in 
the largest merging component, and has been producing stars at a 
rate of 0.8 M yr" 1 . 

Finally, we do not detect the previously reported 'strong' WR 
signatures in UM 448, on a single-spaxel, regional or global scale 
(Masegosa et al. 1991; Schaerer et al. 1999). Using the 4686 A fea- 
ture we tentatively estimate that ~4000 WN stars may be present 
within UM 448. In consideration of the physical size of increased 
N/O seen within UM 448 (~0.6 kpc 2 ) we find that WR stars 
alone cannot be responsible for the region of N-enrichment within 
UM 448. Thus being, another chemical evolution scenario (e.g. in- 
flow of metal-poor gas due to interaction or merger scenario, as the 
morphology and dynamics of UM 448 suggest) is required to ex- 
plain this enhancement. Since there is no strong evidence in support 
of these scenarios (although broad component emission is sugges- 
tive of an outflow, the N-enrichment is seen in the narrow com- 
ponent only), more observations, at higher sensitivity, are needed 
of this intriguing system. In relation to other high N/O BCGs, this 
chemodynamical study supports its membership within the anoma- 
lous subset and, along with other spatially resolved spectroscopic 
studies of high N/O galaxies, demonstrates the benefits that this 
method can have in both isolating the enrichment and aiding our 
understanding of the processes that cause it. 
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APPENDIX A: REGIONAL OBSERVED FLUXES AND 
LINE INTENSITIES 

Tables A1-A3 list the observed fluxes and de-reddened line intensi- 
ties derived from summed spectra over the regions defined in Fig. 2. 
Fluxes are given for each separate velocity component (C1-C3), 
along with its de-reddened line intensity and FWHM. 

This paper has been typeset from a TgX/ ETgX file prepared by the 
author. 
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Table Al. UM 448 regional fluxes and de-reddenened line intensities (both relative to F(H/3)=I(H/3)=1()0) and FWHMs (in km s _1 ) for summed spectra over 
region 1 (as defined in Fig. 2). Line fluxes were extinction-corrected using the c(H/3) values shown at the bottom of the table, calculated from the relative Ha, 
H/3 and H7 fluxes. F(H/3) is in units of X 10~ 14 ergs s^ 1 cm~ 2 arcsec~ 2 . 
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100.26±5.89 


86.57±1.55 


67.05±3.16 


65.31 ±4.59 


191.52±4.09 


[Nil] A6548 


23.57±1.59 


18.37±1.36 


81.05±2.43 


12.88±1.07 


9.09±0.82 


165.02±5.90 


Ha 


370.44± 11.04 


290.03± 12.58 


103.56±0.57 


415.93± 12.45 


295.55± 13.99 


278.17±2.12 


[N 11] A6584 


97.70±2.92 


75.84±3.27 


92.68±0.88 


31.35±1.13 


22.01±U2 


253.20±4.52 


He 1 A6678 


9.41±0.33 


7.23±0.34 


102.59±2.28 








[S 11] A6716 


53.03±2.02 


40.55±1.97 


85.70±1.06 


29.47±1.24 


20.26±1.11 


180.53±2.77 


[S 11] A6731 


43.70±1.59 


33.42±1.58 


85.52±1.06 


18.79±0.91 


12.92±0.77 


180.14±2.76 


c(H/3) 




0.34 ± 0.02 






0.48 ± 0.02 




F(H/3) 




1.95 ± 0.05 






2.72 ± 0.07 





Table A2. Same as for Table Al, for summed spectra over region 2 (as defined in Fig. 2). 
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CI 






C2 






Fa 


Ia 


FHWM 


Fa 


l \ 


FHWM 


[0 II] A3727 


37620±7.41 


506.95±52.62 


199.37±0.73 








[0 II] A3729 


487.54±9.55 


656.99±68.17 


212.46±15.14 








[Nem] A3868 


65.57±4.66 


85.82±10.52 


213.38±19.76 








H8+He I A3888 


55 24±5.09 


72.05±9.78 


221.94±2.30 








HS 


9.51 ±2. 11 


11.76±2.85 


95.75±9.88 


26.11±1.09 


32.67±1.87 


157.59±5.69 




.1 1 . o ) It J . J J 


Af. 9^-1-^ 7T 
H-U.ZJ It J. / J 


J 1 D.L JZt — O.Uo 






TJfi lf,-\-00 96 
ZZU. /O^ — l. — U 


[0 III] A4363 


1.45±0.84 


1.67±0.97 


74.34±37.30 


5.33±0.07 


6.19±0.24 


455.65±67.72 


He I A4471 


13.64±1.33 


15.25±2.01 


226.35±22.75 








H/3 


100.00±2.71 


100.00±8.56 


118.95±8.93 


100.00±0.96 


100.00±3.46 


233.89±7.09 


[0 III] A4959 


97.79±2.17 


95.09±7.82 


94.79±3.52 


89.70±0.97 


87.09±2.98 


255.45±4.99 


[N II] A6548 


60.50±2.56 


41.76±2.91 


188.97±2.54 








Ha 


414.10±8.58 


287.88±17.14 


111.69±4.90 


417.95±3.94 


284.58±7.03 


249.28±5.71 


[N II] A6584 


132.16±2.76 


90.71 ±5.33 


162.09±7.08 


35.66±3.51 


23.95±2.42 


333.30±17.45 


He I A6678 


14.28±0.71 


9.66±0.71 


200.80±7.71 








[S II] A6716 


176.68±4.77 


1 18.61 ±7. 10 


193.47±0.69 








[S II] A6731 


135.17±4.06 


90.75±5.56 


193.01 ±0.69 








c(H/3) 




0.51 ± 0.04 






0.53 ± 0.01 




F(H/3) 




1.67 ±0.03 






3.54 ± 0.02 





C3 




Fa 


Ia 


FHWM 


[On] A3727 








[O II] A3729 








[Nem] A3868 








H8+Hel A3888 








HS 








H 7 


41.93±3.60 


47.58±5.45 


106.30±20.01 


[O III] A4363 








He I A4471 








H/3 


100.00±5.80 


100.00±8.88 


84.49±6.54 


[O III] A4959 


129.25±8.12 


126.20± 11.47 


86.50±4.99 


[N II] A6548 








Ha 


401.60±16.70 


294.40± 18.32 


86.10±4.71 


[N II] A6584 


77.15±3.50 


55.94±3.60 


91.24±6.78 


Hel A6678 








[S II] A6716 








[S II] A6731 








cW) 




0.43 ± 0.03 




F(H/3) 




0.54 ± 0.02 





Table A3. Same as for Table Al, for summed spectra over region 3 (as defined in Fig. 2). 
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ABSTRACT 

Using VLT/FLAMES optical integral field unit observations, we present a detailed study of 

r ~~^ UM 448, a nearby Blue Compact Galaxy (BCG) previously reported to have an anomalously 

high N/O abundance ratio. NTT/SuSI2 images reveal a morphology suggestive of a merger of 

\^ two systems of contrasting colour, whilst our Ha emission maps resolve UM 448 into three 

• separate regions that do not coincide with the stellar continuum peaks. UM 448 exhibits com- 

'-^ plex emission line profiles, with most lines consisting of a narrow (FWHM<100 kms -1 ), 

central component, an underlying broad component (FWHM^ 150-300 kms -1 ) and a third, 
narrow blue-shifted component. Radial velocity maps of all three components show signs of 
S-j solid body rotation across UM 448, with a projected rotation axis that correlates with the con- 

tinuum morphology of the galaxy. A spatially-resolved, chemodynamical analysis based on 
the [O in] AA4363, 4959, [N n] A6584, [S n]AA6716, 6731 and [Ne in] A3868 line maps, is 
presented. Whilst the eastern tail of UM 448 has electron temperatures (T ) that are typical 

y— \ of BCGs, we find a region within the main body of the galaxy where the narrow and broad 

J> [O III] A4363 line components trace temperatures differing by 5000 K and oxygen abundances 

t-H differing by 0.4 dex. We measure spatially resolved and integrated ionic and elemental abun- 

dances for O, N, S and Ne throughout UM 448, and find they do not agree, possibly due the 
flux- weighting of T c from the integrated spectrum. This has significant implications for abun- 
dances derived from long-slit and integrated spectra of star-forming galaxies in the nearby and 

Q\ distant universe. A region of enhanced N/O ratio is indeed found, extended over a ^0.6 kpc 2 

region within the main body of the galaxy. Contrary to previous studies, however, we do not 

04 find evidence for a large Wolf-Rayet population, and conclude that WR stars alone cannot be 

responsible for producing the observed N/O excess. Instead, the location and disturbed mor- 
phology of the N-enriched region suggests that interaction-induced inflow of metal-poor gas 

• may be responsible. 

x 

5_i Key words: galaxies:abundances - galaxies: individual (UM 448) - galaxies: dwarf - galax- 

ies:kinematics and dynamics - galaxies interactions - stars: Wolf-Rayet 



1 INTRODUCTION 

Observations of Blue Compact Galaxies (BCGs) in the nearby Uni- 
verse provide a means for studying chemical evolution and star 
formation (SF) processes in chemically un-evolved environments. 
Their low metallicities make them attractive as analogues to the 
young 'building-block' galaxies thought to exist in the high-z pri- 
mordial universe (Searle & Sargent 1972). BCGs are thought to 
have experienced only very low-level star formation in the past, 
but are currently undergoing bursts of SF in relatively pristine 
environments, ranging from 0.02-0.5 solar metallicity (Kunth & 
Ostlin 2000). Thus understanding the chemical evolution within 

* E-mail: bjames@stsci.edu 



BCGs can impact upon our understanding of primordial galaxies 
and galaxy evolution in general. 

The distribution of nitrogen as a function of metallicity in SF 
galaxies has often been a cause for debate, especially at interme- 
diate metallicities (7.6< 12+log(0/H) <8.3). Within this metallic- 
ity range the N/O abundance ratio indicates that nitrogen behaves 
neither as a primary element (i.e. N/O independent of O/H) nor 
as a secondary element (i.e. N/O proportional to O/H); instead a 
large scatter is observed. In particular, several BCGs exist whose 
N/O abundance ratios are up to three times those expected for their 
O/H values, as highlighted by e.g. Pustilnik et al. (2004). Several of 
these galaxies have been the focus of our spectroscopic studies us- 
ing integral field unit (IFU) observations (James et al. 2009, 2010), 
aimed at understanding the source of such abundance peculiarities. 
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James et al. (2009) performed a detailed abundance analysis 
of the N/O-outlier BCG, Mrk 996, whose X-ray properties have 
been studied by Georgakakis et al. (2011). Previous long-slit ob- 
servations of this galaxy (Thuan et al. 1996; Pustilnik et al. 2004), 
had shown it to have N/O ratio values at least ~4 times higher 
than those of other BCGs of similar metallicity. Using IFU obser- 
vations, James et al. (2009) performed a spatially resolved multi- 
component emission line analysis of Mrk 996, which revealed a 
twenty-fold N/O and N/H enrichment factor in the broad compo- 
nent (400-500 km/s full width at half maximum - FWHM) versus 
the narrow component line emission. The narrow component gas 
on the other hand has a N/O ratio which is typical for the galaxy's 
metallicity. In this respect, Mrk 996 is similar to NGC 5253, a star- 
burst galaxy for which an enhanced N/O ratio has also been iso- 
lated to broad component gas (Lopez-Sanchez & Esteban 2010). 
In Mrk 996, the spatial morphology of the broad line component 
is highly correlated with Wolf-Rayet (WR) wind-feature emission 
arising from a population of WNL stars, strongly suggesting the 
enhanced nitrogen to be the result of chemical enrichment by the 
N-rich winds of WN-type WR stars. This supports the postulate of 
nitrogen self-enrichment due to WR stars that had been suggested 
by Walsh & Roy (1987) and Pagel et al. (1986) to explain elevated 
N/O levels. 

James et al. (2010) also presented an integral field spectro- 
scopic study of the interacting dwarf galaxy UM 420 for which 
the long-slit spectroscopic study of Izotov & Thuan (1999) had re- 
ported log(N/0) = — 1 .08. By performing a spatially resolved anal- 
ysis of its physical and chemical conditions, James et al. (2010) 
found an average ratio applicable to the whole galaxy of log(N/0) 
= —1.64, i.e. 3.6 times lower than previous estimates. The VLT 
observations of James et al. (2010) did not confirm previous detec- 
tions of WR spectral features in UM 420, highlighting the difficul- 
ties involved in identifying WR stellar populations in distant dwarf 
galaxies. 

In this paper, we present IFU observations of another mem- 
ber of this subset of potentially N/O-anomalous galaxies, UM 448 
(UGC 06665, Mrk 1304). With a reported factor of ~3 nitrogen 
excess (Bergvall & Ostlin 2002), this galaxy presents an excel- 
lent opportunity to further our understanding of peculiar N/O ra- 
tios in some BCGs and the mechanisms responsible. UM 448 is 
a complex system, with a disturbed morphology possibly due to a 
merger event (Dopita et al. 2002). Thus the need for spatially re- 
solved spectroscopy is two-fold: firstly in determining the range of 
physical properties and chemical compositions that are present, and 
secondly, in spatially relating any chemical anomalies to environ- 
mental factors (e.g. WR populations). 

UM 448 is a starburst galaxy (Keel & van Soest 1992; Kew- 
ley et al. 2001) consisting of various star forming knots, diffuse 
Ha emission in the north, and a SE tidal extension (Dopita et al. 
(2002)). Broad He II A4686 emission was detected by Masegosa 
et al. (1991) and subsequently by Guseva et al. (2000), who re- 
ported an atypical weakness in this line in comparison with other 
WR-signature lines. As with Mrk 996 and NGC 5253, a WR popu- 
lation may well be the source of N-enhancement within this BCG. 

In this paper we present integral field spectroscopy (IFS) ob- 
servations obtained with the VLT-FLAMES IFU. Crucially, this in- 
strument is sufficiently blue-sensitive to detect the [O II] A3727 
doublet line, meaning that we are not required to adopt uncer- 
tain ionisation correction factors to determine oxygen abundances. 
These data afford us a new spatiokinematic 3-D view of UM 448, 
with the spatial and spectral resolution to undertake a full multi- 
velocity component analysis. We adopt a distance of 75.9 Mpc for 



UM 448 from the velocity measured here of +5579±5 km s 
(z = 0.018660 ± 1.6 x 10~ 5 ), with a Hubble constant of 
H =73.5 km s" 1 Mpc" 1 (DeBernardis et al. 2008). 



2 OBSERVATIONS AND DATA REDUCTION 

Observations of UM 448 were obtained with the Fibre Large Ar- 
ray Multi Element Spectrograph, FLAMES (Pasquini et al. 2002) 
at Kueyen, Telescope Unit 2 of the 8.2 m VLT at ESO's observatory 
on Paranal, in service-mode on dates specified in Table 1. Obser- 
vations were made with the Argus Integral Field Unit (IFU), with 
a field-of-view (FoV) of ll'.'5x7'.'3 and a sampling of O'.'52/lens. 
In addition to science fibres, Argus has 15 sky-dedicated fibres that 
simultaneously observe the sky and that were arranged to surround 
the IFU FoV at distances of ~35-58". The positioning of the IFU 
aperture for UM 448 is shown overlaid on SuSI2 images 1 in Fig. 1. 

Four different low-resolution (LR) gratings were utilised: 
LR1 (AA3620-4081, 24.9±0.3 kms" 1 FWHM resolution), 
LR2 (AA3964-4567, 24.7±0.1 kms" 1 ), LR3 (AA4501-5078, 
24.6±0.3 kms" 1 ) and LR6 (AA6438-7184, 21.9±0.4 kms" 1 ). 
This enabled us to cover all the important emission lines needed for 
an optical abundance analysis, as detailed in Section 3. The photo- 
metric conditions, airmass and exposure time for each data set are 
detailed in Table 1. In addition to the science frames, spectropho- 
tometric standard star observations (Table 1), as well as continuum 
and ThAr arc lamp exposures were obtained. 

The data were reduced using the girBLDRS pipeline (Blecha 
& Simond 2004) 2 in a process which included bias removal, lo- 
calisation of fibres on the flat-field exposures, extraction of indi- 
vidual fibres, wavelength calibration and rebinning of the ThAr 
lamp exposures, and the full processing of science frames which 
resulted in flat-fielded, wavelength-rebinned spectra. The pipeline 
recipes 'biasMast', 'locMast', 'wcalMast', and 'extract' were used 
(cf. Tsamis et al. 2008). The frames were then averaged using the 
'imcombine' task of IRAF 3 which also performed the cosmic ray 
rejection. The flux calibration was performed within IRAF using 
the tasks CALIBRATE and STANDARD. Spectra of the spectropho- 
tometric stars quoted in Table 1 were individually extracted with 
girBLDRS and the spaxels containing the stellar emission were 
summed up to form a ID spectrum. The sensitivity function was 
determined using SENSFUNC and this was subsequently applied 
to the combined UM 448 science exposures. The sky subtraction 
was performed by averaging the spectra recorded by the sky fi- 
bres and subtracting this spectrum from that of each spaxel in the 
IFU. Custom-made scripts were then used to convert the row by 
row stacked, processed CCD spectra to data cubes. This resulted in 
one science cube per grating, i.e. four science cubes for the source 
galaxy. 

1 SuSI2 images were obtained from the ESO archive, with total exposure 
times of 1200, 1200, 1800 and 2700 s, average airmasses of 1.34, 1.28, 1.22 
and 1.17 and image qualities of 0.75", 0.93", 0.87" and 0.56" FWHM for 
the R-, V-, B- and l/-band images, respectively. The average image qual- 
ity of the images is ~0.78" FWHM. WCS information was assigned via 
position matching with GSC2. Other than aligning and median-combining 
the separate frames using IRAF's IMALIGN and IMCOMBINE, respectively, 
no additional reduction was performed on the images. 

2 The girBLDRS pipeline is provided by the Geneva Observatory 
http://girbldrs.sourceforge.net 

3 IRAF is distributed by the National Optical Astronomy Observatory, 
which is operated by the Association of Universities for Research in As- 
tronomy 
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Figure 1. SuSI2 (Superb-Seeing Imager, installed on the ESO 3.5m New Technology Telescope (NTT) at La Silla Observatory, Chile) colour-composite 
images of UM 448 (red, i?-band; green, V-band; blue B-band), with the 11,5" X7.3" FLAMES IFU aperture overlaid. Each image has been scaled to 
enhance different features within UM 448. The left-panel shows two large knots of continuum or ionised gas emission that extend away from the galaxy in a 
south-westerly direction, outside the IFU FoV. The right-panel shows the two distinct portions of the galaxy, a blue eastern tail and a reddish knot of stellar 
continuum. SuSI2 images were obtained from the ESO archive, observed as part of programme ()78.B-0358(A). 



Table 1. FLAMES IFU observing log 

Date Grism Wavelength Range Exp. time Avg. Airmass FWHM seeing Standard Star 







(A) 


(s) 




(arcsec) 




24/05/2009 


LR1 


3620-4081 


3x233 


1.11 


0.85 


LTT7379 




LR2 


3964-4567 


3x300 


1.12 


0.87 


LTT7379 




LR3 


4501-5078 


3x230 


2.24 


0.80 


LTT7379 


18/06/2009 


LR1 


3620-4081 


3x233 


1.12 


0.47 


Feigel 10 




LR2 


3964-4567 


3x300 


1.13 


0.42 


Feigel 10 




LR3 


4501-5078 


3x230 


1.16 


0.50 


Feigel 10 


28/07/2009 


LR6 


6438-7184 


3x400 


1.75 


0.47 


EG21 


26/01/2010 


LR1 


3620-4081 


3x233 


1.15 


0.54 


LTT4364 




LR2 


3964-4567 


3x300 


1.13 


0.53 


LTT4364 




LR3 


4501-5078 


3x230 


1.11 


0.63 


LTT4364 


16/03/2010 


LR1 


3620-4081 


3x233 


1.10 


0.57 


LTT4364 




LR2 


3964-4567 


3x300 


1.11 


0.50 


LTT4364 




LR3 


4501-5078 


3x230 


1.12 


0.56 


LTT4364 


08/07/2010 


LR6 


6438-7184 


8x600 


1.63 


0.77 


Feigel 10 


09/07/2010 


LR1 


3620-4081 


3x233 


1.50 


0.77 


Feigel 10 




LR2 


3964-4567 


3x300 


1.60 


0.72 


Feigel 10 




LR3 


4501-5078 


3x230 


1.74 


0.66 


Feigel 10 



Observations of an object's spectrum through the Earth's at- 
mosphere are subject to refraction as a function of wavelength, 
known as differential atmospheric refraction (DAR). The direction 
of DAR is along the parallactic angle at which the observation is 
made. Following the method described in James et al. (2009), each 
reduced data cube was corrected for this effect using the algorithm 
created by Walsh & Roy (1990) ; this proceedure calculates the 
fractional pixel shifts for each monochromatic slice of the cube rel- 
ative to a fiducial wavelength (i.e. a strong emission line), shifts 
each slice with respect to the orientation of the slit on the sky and 
the parallactic angle and re-combines the DAR-corrected cube. All 
pixel shifts were < 1 pixel in size, indicating uncertainties due to 
extinction or temperature gradients are negligible. 



3 MAPPING OF LINE FLUXES AND KINEMATICS 



The distribution of H I Balmer line emission across UM 448 which 
is indicative of current massive star formation activity, and the mor- 
phology of the system as it appears on the SuSI2 images, were used 
to define three regions of roughly similar Balmer line luminosity. 
Our analysis of the galaxy is based on properties averaged over 
each of these three areas whose boundaries are displayed on the 
Ha map in Figure 2. Regions 1 and 2 are both contained in the blue 
eastern tail of UM 448 as seen on Fig. 1 (right), containing the bulk 
of Balmer line emission. Region 3 encompasses the interface of the 
blue eastern tail with the redder western part of the galaxy, and fur- 
ther contains the peak of the continuum emission in the vicinity of 
Ha (see Fig. 12). Whilst each of these regions shows a relatively 
small variation in surface brightness, the composite structure of the 
emission lines can vary considerably as illustrated by the Ha emis- 
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Figure 2. Ha emission per 0'.'52x(X'52 arcsec 2 spaxel. Three regions were defined across UM 448 over which summed spectra were extracted for analysis. 
On the right the Ho line profile from each summed region is shown, together with a two-component Gaussian fit. Blue dashed lines demarcate the wavelength 
range used for summing the flux. North is up and east is to the left. 





Fa 




[On] A3727 


98.02 ± 1.18 


125.40 ± 6.20 


[Oil] A3729 


135.65 ± 1.43 


173.55 ± 8.52 


[Nelll] A3868 


23.41 ± 1.71 


29.24 ± 2.54 


H8+Hel A3888 


15.47 ± 2.24 


19.26 ± 2.93 


HS 


21.57 ± 1.21 


25.70 ± 1.85 


H7 


43.33 ± 0.83 


48.96 ± 2.31 


[OIII] A4363 


1.90 ± 0.19 


2.14 ± 0.23 


He I A4471 


4.50 ± 0.49 


4.94 ± 0.57 


H/3 


100.00 ± 0.99 


100.00 ± 3.94 


[Olll] A4959 


91.86 ± 1.66 


89.76 ± 3.71 


[Nil] A6548 


15.83 ± 0.52 


11.66 ±0.49 


Ha 


407.05 ± 13.27 


301.47 ± 12.62 


[N II] A6584 


58.55 ± 0.51 


42.91 ± 1.17 


He 1 A6678 


2.86 ± 0.06 


2.07 ± 0.07 


[Sll] A6716 


49.17 ± 0.44 


35.38 ± 0.94 


[S ll] A6731 


36.07 ± 0.34 


25.95 ± 0.70 


c(H/3) 


0.42 ± 0.02 


F(H/3) 


24.48 ±0.17 



Table 2. UM 448 fluxes and de-reddenened line intensities (both relative 
to F(H/3)=I(H/3)=100) for summed spectra over all of UM 448. Line fluxes 
were extinction-corrected using the c(H/3) values shown at the bottom of 
the table, calculated from the relative Ha, H/3 and H7 fluxes. The line in- 
tensities listed here were used for summed-spectra T c and iV c diagnostics 
and regional ionic abundance calculations listed in Table 3. 



nents in the majority of emission lines seen in the 300 spectra of 
UM 448 across the IFU aperture. Following the methods outlined in 
James et al. (2009) and Westmoquette et al. (2007), we utilised an 
automated fitting proceedure called PAN (Peak ANalysis; Dimeo 
2005), that allowed us to read in multiple spectra simultaneously, 
interactively specify initial parameters for a spectral line fit and 
then sequentially process each spectrum, fitting Gaussian profiles 
accordingly. The output consists of the \ 2 value for the fit, the fit 
results for the continuum and each line's flux, centroid and FWHM, 
and the respective errors of these quantities which were propagated 
in all physical properties derived afterwards 4 . 

Even though some low S/N lines may in reality be composed 
of multiple emission components, it is often not statistically signif- 
icant to fit anything more than a single Gaussian. In order to rig- 
orously determine the optimum number of Gaussians required to 
fit each observed profile, we followed the likelihood ratio method 
outlined by Westmoquette et al. (201 1, their Appendix A). S/N ra- 
tio maps were first made for each emission line, taking the ratio 
of the integrated intensity of the line to that of the error-array pro- 
duced by the FLAMES pipeline, on a spaxel-by-spaxel basis. Sin- 
gle or multiple Gaussian profiles were then fitted to the emission 
lines, restricting the minimum FWHM to be the instrumental width 



sion line profiles displayed alongside the monochromatic line map 
shown in Figure 2. 

The high spectral resolution and signal-to-noise (S/N) ratio 
of the data allow the identification of multiple velocity compo- 



4 The uncertainties output by the pipeline were used as the starting point of 
the enor analysis; they also include uncertainties from the flat-fielding and 
the geometric correction of the spectra and not just photon-noise (Blecha & 
Simond 2004). 



(see Section 2). Suitable wavelength limits were defined for each 
line and continuum level fit. Further constraints were applied when 
fitting the [S II] and [O II] doublets: the wavelength difference 
between each line in the doublet was taken to be equal to the red- 
shifted laboratory value when fitting the velocity component, and 
their FWHMs were set equal to one another. 

The majority of high S/N ratio line profiles in the spectra 
of UM 448 are optimally fitted with a single narrow Gaussian 
(FWHM<100 kms , component CI hereafter), a single under- 
lying broad Gaussian component (FWHM 150-300 kms - , C2), 
and a third narrow Gaussian which can appear either red- or blue- 
shifted with respect to CI (C3 hereafter). In consideration of the 
different 'features' of each component (i.e. FWHM and/or veloc- 
ity), we consider each component as being real, i.e. arising from 
gas with different physical conditions. Towards the outer regions 
of the galaxy, where the S/N ratio decreases, the lines are optimally 
fitted with only CI and C2 velocity components. Low S/N ratio 
emission lines (e.g. H8+He I A3888) were optimally fitted with a 
single narrow Gaussian profile (CI). Figure 3 shows the complex 
profile of the Ha line, along with the spatial distribution of its con- 
stituent velocity components. 

Tables A1-A3 list measured FWHMs and observed and de- 
reddened fluxes for the fitted Gaussian components of the detected 
emission lines across the three regions of massive star formation 
defined in Fig. 2. The fluxes are for spectra summed over each re- 
gion and are quoted relative to the flux of the corresponding H/3 
component. They were corrected for reddening using the galactic 
reddening law of Howarth (1983) using c(H/3) values derived from 
the Ha/H/3 and H7/H/3 line ratios of their corresponding compo- 
nents, weighted in a 3:1 ratio °, respectively, in conjunction with 
theoretical Case B ratios from Storey & Hummer (1995). An ex- 
tinction map was also derived using the the Ha/H/3 and H7/H/3 
emission line ratio maps and is shown in Fig. 4. This was com- 
puted using the integrated line profiles (i.e. C1+C2+C3) to obtain 
a global extinction view of the galaxy. 

A Milky Way reddening of E(B - V) = 0.025 mag in the di- 
rection to UM 448 is indicated by the extinction maps of Schlegel 
et al. (1998), corresponding to c(H^)=0.036. Total c(H^) values ap- 
plicable to the galaxy and its individual emission line components 
in each main aperture-region are listed in their respective tables 
(Tables A1-A3). 

3.1 Ha maps 

The Ha surface brightness of UM 448 (Fig. 3a) displays a hor- 
izontal 'S'-shaped morphology, extending a total of ~ 10" east 
to west, almost symmetrically about the FoV center. The lower 
arm of the galaxy extends ~ 5" to the east and contains the 
galaxy's peak Ha emission (corresponding to designated region 
1). Based on our fitting method, the average surface brightness 
within region 1 can be de-composed to ~ 1.8 x 10" 14 and ~ 
2.8 x 10 -14 ergs cm -2 s -1 arcsec -2 for line components CI and 
C2 respectively. The upper arm is seen to extend ~ 5" in the oppo- 
site direction and displays an overall slightly lower surface bright- 
ness of ~ 1.2 x 10" 14 ergs cm" 2 s -1 arcsec -2 for both compo- 
nents. Component C3 appears in three separate regions; two regions 

5 An exception was made for region 3 where the low S/N ratio of H7 pro- 
hibited the robust decomposition into the three components observed in Ha 
and H/3. Therefore, only the Ha/H/3 line ratio was used for the calculation 
of c(H/3) for this region. 



A 3D Chemodynamical Study of UM 448 5 




-4 -2 2 4 
arcsecs 



0.00 0.19 0.38 0.56 0.75 0.94 
E(B-V) 

Figure 4. Map of E(B-V) for UM 448, derived from the integrated flux of 
H I emission lines. Integrated light Ho emission line contours are overlaid. 
North is up and east is to the left. 



either side of the central boundary region, i.e. between the two stel- 
lar components (see Fig. 1), and a third region along the most west- 
ern part of the upper arm. C3 had an average surface brightness of 
~ 4.4 x 10 -14 ergs cm" 2 s" 1 arcsec" 2 . The total integrated Ha 
flux for UM 448 is 9.24x 10 -13 ergs cm -2 s -1 , which is in good 
agreement with the value of ~ 11.5 x 10 -13 ergs cm -2 s -1 de- 
rived by Dopita et al. (2002) from narrow band Ha imaging. Rela- 
tive line intensities are listed in Tables A1-A3 for summed spectra 
over each of the three designated regions, for each separate velocity 
component. 

The Ha velocity map (Fig. 3b) displays a rather similar ra- 
dial velocity structure for both CI and C2 components. Both show 
a smooth velocity distribution between +50 and —180 kms -1 , 
which suggests solid body rotation within the galaxy. These maps 
can be used to derive an axis of rotation for the gaseous component 
of UM448 along position angle, PA ~ 36° (shown as a dashed line 
in Fig. 3b). Interestingly, the alignment of this axis also divides the 
blue eastern 'tail' from the red, western part of the galaxy (Fig. 1). 
The radial velocity structure of component C3 is less straightfor- 
ward than its counterparts. In the two regions near the center of 
UM 448, the velocities range from -100 to +100 kms -1 , whilst 
in the extreme west of the galaxy they are consistently negative 
down to —300 kms -1 . The velocity structure there could be the 
signature of a discrete blue-shifted outflow, which is perhaps re- 
lated to the detached knots of emission present in the SuSI2 images 
~ 10" from the southwest corner of the Argus aperture (Fig. 1) 
- i.e. we are detecting fast-moving gas, kinematically decoupled 
from the main body of UM448. 

The maps of the Ha FWHM show a slight spatial correlation 
with the flux distribution (Fig. 3c). The eastern regions displays 
decreased FWHMs for both the CI and C2 line components, in 
comparison to the upper arm structure: the average FWHM for the 
lower arm is -40-60 kms" 1 for CI and 120-160 kms" 1 for C2, 
which increases to ~90-100 km s" 1 and 200-240 km s" 1 , respec- 
tively, in the upper arm. Interestingly, the region of low FWHM 
is co-spatial with the main peak in flux within the eastern arm. If 
the overall appearance of UM448 is indicative of a close encounter 
or merger of two systems, then it would make sense for the gas 
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Figure 3. Maps of UM 448 in the Ho velocity components C1-C3: the top row shows an example 3-component fit to the Ha line and the structure of the 
velocity component shown in red is displayed below; a) flux in units of 10 — 16 ergs s~ 1 cm -2 arcsec~ 2 per 0'.'52x0'. / 52 spaxel; b) radial velocity relative 
to the heliocentric systemic velocity of +5579 kms —1 (the dashed line refers to the PA~36° rotation axis discussed in the text; c) FWHM corrected for the 
instrumental PSF. See text for details. Maps are shown for S/N > 10 in the integrated line profile. North is up and east is to the left. 



turbulence, and hence the FWHM, to increase as one approaches 
the transition region between the two colliding components (the 
blue/red interface in Fig. 1 right). On the other hand, component 
C3 shows a typical FWHM of 70 km s _1 suggesting an association 
with more orderly gas volumes. 

3.2 Forbidden-line maps 

Emission line maps of UM 448 in the light of the [O III] [N II] 
and [Ne III] lines for the various velocity components are shown 
in Fig. 5. The morphology of [O III] A4959 is similar to that of 
Ha, with the surface brightness of the lower arm being higher by a 
factor of ~2 than that of the upper arm (for component CI). In con- 
trast, [N II] displays a rather constant surface brightness throughout 
either arm, with a peak in emission within our designated region 2. 
For both [O III] and [N II] line component C2 displays a simi- 
lar morphology to the CI counterparts, as was the case with Ha. 



Component C3 of [O III] A4959 has a wider spatial extent than 
its counterpart in Ha. For [N II] component C3 is seen to overlap 
roughly with the central and eastern regions observed for Ha-C3. 
In the central regions C3 shows symmetry approximately about the 
rotation axis defined in Fig 3b; this may a signature of an outflow, 
or it simply reflects the highly disturbed kinematics in the "collision 
zone" of the merging galactic components at the Region l/Region2 
aperture interface. [O III] A4363 has also been detected in veloc- 
ity components CI and C2 whose peak surface brightness is found 
within region 3 (see Fig. 10 for an example fit to the summed spec- 
tra across this region), whereas A4959 peaks in region 1. This is 
mirrored in the high electron temperature of region 3 versus region 
1 (see below). The [Ne III] A3868 emission also peaks within re- 
gion 1, with more diffuse emission extending throughout regions 2 
and 3. 

UM 448 is shown in the light of the [O II] and [S II] dou- 
blet lines in Fig. 6. Due to their relatively low S/N ratio (compared 
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Figure 5. Dereddened intensity maps of UM 448 in [O III] A4959, [N II] A6584, [O III] A4363 and [Ne III] A3868 for the separate velocity components CI, 
C2 and C3. [O III] A4363 is shown after being re-binned by 1,5x1.5 spaxels and re-mapped onto the original grid (see the text for details). Only spaxels with 
S/N> 10 are shown for [O III] A4959 and [N II] A6584, >5 for [Ne III] A3868 and > 1 for [O III] A4363. North is up and east is to the left. 



to the H I Balmer lines and [O III] A4959), both doublets could be 
decomposed to only the CI and C2 velocity components (Fig. 7). 
The [O II] doublet, with a rest-wavelength separation of 2.783 A 
was particularly hard to decompose into a narrow+broad compo- 
nent since the large FWHM of the component C2 would effectively 
merge the doublet into a single line. To alleviate this, each spaxel 
with a broad component [O II] line detection was examined and 
only selected when there was sufficient S/N ratio to allow for a de- 
blended fit (see Fig. 7a). The larger wavelength separation of the 
[S II] doublet allowed the broad component (C2) to be also fit- 
ted in relatively low S/N ratio spaxels (e.g. Fig. 7b), and was thus 
mapped more extensively than [O II] (compare the lower panels 
of Fig. 6a, b). The [O II] A3727+A3729 emission has a rather dis- 
jointed morphology compared to the integrated Ha surface bright- 
ness. The ratio of the [O II] doublet components commonly used 
as a density diagnostic, (Fig. 6, left panels) displays negligible 
structure in relation to the overlaid Ha. In contrast to [O II], the 
[S II] A6716+A6731 doublet has a morphology that is more akin 
to Ha with three distinct peaks in flux. The A6716/A6731 ratio 
(Fig. 6, right panels), which is also used as a density diagnostic, 
decreases in regions that align with the peaks in Ha emission, and 
this is indicative of increased electron density there. 



4 ELECTRON TEMPERATURE AND DENSITY 

The [O III] (A5007 + A4959)/A4363 intensity ratios were used to 
determine electron temperatures. Two emission line doublets are 
available to determine the electron density, N a ; namely the [O II] 
A3726/A3729 and [S II] A6716/A6731 ratios. T c and N c values 
were computed using IRAF'S ZONES task in the NEBULA pack- 
age. This makes use of technique that derives both the tempera- 
ture and the density by making simultaneous use of temperature- 
and density-sensitive line ratios. Atomic transition probabilities for 
2+ andO+ were taken from Wiese et al. (1996) whilst collision 
strengths were taken from Lennon & Burke (1994) and McLaugh- 
lin & Bell (1993), respectively . For [S II] we used the transitional 
probabilities of Keenan et al. (1993) and the collisional strengths of 
Ramsbottom et al. (1996). 

Where possible T e and N e values were derived for each ve- 
locity component within each of the aperture regions defined in 
Fig. 2 and are listed in Table 3. The mean electron densities from 
the [S II] and [O II] ratios agree to within ~20%; the [O II] density 
map is however significantly more noisy. The [S II] densities were 
therefore adopted for the analysis. 

T e and N a maps computed for line components CI and C2 are 
shown in Fig. 8, and mean values for each of the aperture regions 
1-3 of Fig. 2 are listed in Table 3. Fig. 8 (left) shows evidence for 
a small amount of T c variation throughout UM 448. Aperture re- 
gions 1 and 2 show temperatures of ~ 12,000-13,000 K. Region 3 
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Figure 6. Dereddened emission line maps of UM 448 in doublet ratioes [O II] A3727/A3729 (top-panel) and [O II] A6716/A6731(bottom-panel), in the separate 
Cl(narrow) and C2 (broad) velocity components, with integrated Ha emission line flux overlaid in contours. North is up and east is to the left. 



shows temperatures of 15,300 K for the narrow line component 
CI and 10,200 K for the broad line component C2. iV B variations 
are also apparent and denser areas are correlated with peaks in Ha 
emission. Regions 1 and 2 have an average 7Y C ~120-150 cm -3 , 
whereas region 3 is only marginally denser with 200 cm -3 . 



5 CHEMICAL ABUNDANCES 

Ionic abundance maps relative to H + were created for the N + , + , 
2+ , Ne 2+ and S + ions, using the AA6584, 3727+3729, 4959, 
3868 and 67 17+673 1 lines respectively. Examples of the flux maps 
used in the derivation of these ionic abundance maps can be seen 
in Figure 5. Abundances were calculated using the ZONES task 
in IRAF, using the respective T and 7V C maps described above, 
with each FLAMES spaxel treated as a distinct 'nebular zone' with 
its own set of physical conditions. These results were compared 
with those obtained using the multilevel ionic equilibrium program 
EQUIB (originally written by I. D. Howarth and S. Adams at Univ. 
College London), which comprises a different set of atomic data 
tables. Typical differences in the abundance ratios were found to be 
5% or less. To check the validity of using T c ([0 III]) for the singu- 
larly ionized elements (N + , S + , and + ), we also computed abun- 
dance maps useing a T c ([0 II]) map (created using the relationship 
with T c ([0 III]) described in Izotov et al. (2006)). Discrepancies 
between the two temperatures were found to be low (~ 100-400 K) 
and typical differences in ionic abundance were ~2-7%, These 



errors are lower than those caused by uncertainties in T c ([0 III]) 
alone and we therefore do not include it as a source of error. 

Ionic nitrogen, neon and sulphur abundances were converted 
into N/H, Ne/H and S/H abundances using ionisation correction 
factors (ICFs) from Kingsburgh & Barlow (1994).The O/H abun- 
dance was obtained by adding the 2+ /H + and + /H + abundance 
maps (i.e. assuming that the O /H + ratio is negligible). Since the 
[S III] A6312 line was not within the wavelength range of our data, 
the S 2+ /H + abundance was estimated using the empirical relation- 
ship between the S 2+ and S + ionic fractions from the corrected 
equation A38 of Kingsburgh & Barlow (1994), namely S 2+ /S + 
= 4.677 x (0 2+ /0+) 433 . The S/H results should therefore be 
considered with caution, but they are not at variance with those 
obtained using the ICF prescriptions of Izotov et al. (2006). O/H, 
N/H and log(N/0) abundance maps are shown in Fig. 9 and are 
described below. Regional abundances derived directly from tak- 
ing error-weighted averages over the abundance maps are listed in 
Table 3. 

5.1 Maps 

The elemental oxygen abundance map shown in the middle panel of 
Figure 9 displays a single peak whose location correlates with the 
peak in Ha relating to region 1 . Minimal abundance variations are 
seen across UM 448, with a slightly decreasing oxygen abundance 
as one moves from regions 1-3, with an overall average abundance 
for CI of 12 + log(0/H) = 8.15 ± 0.10. Adopting a solar oxy- 
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Figure 8. [O III] electron temperature maps of UM 448 for CI (narrow) and C2 (broad) velocity components. Overlaid circles represent the size of the T a 
uncertainty within each spaxel (percentage errors are scaled by the 0.52" diameter of a spaxel). North is up and east is to the left. 
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Figure 7. Example single-spaxel spectra and fits to the 
[O II] AA3727,3729 and [S II] AA67 16,6731 lines. Flux is in units of 
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gen abundance of 8.7 1±0. 1 relative to hydrogen (Scott et al. 2009), 
this corresponds to ~0.28 Zq. This is comparable to the value ob- 
tained by Izotov & Thuan (1998, IT98 hereafter), who measured 
an abundance of 12 + log(0/H) = 7.99±0.04. The C2 veloc- 
ity component displays a higher oxygen abundance than the Cl- 



component gas within region 3 by ~0.4 dex, which, judging from 
the agreement between CI and C2 [O III] 4959 A flux (Fig. 5) is 
mostly attributable to the T c being ~5000 K cooler (if the CI re- 
gion was instead much denser but of similar oxygen abundance to 
the gas sampled by the C2 line component one might expect the 
A4959/H/3 CI ratio to be much smaller something which is not ob- 
served - Table A3). 

The upper panel of Fig. 9 shows the N/H abundance across 
UM 448. The CI peak in nitrogen abundance is offset from that of 
oxygen, peaking at the location of region 2 with an average N/H 
ratio of 1.04±0.01 x 10~ 5 . For region 1 it is slightly lower at 
9.40±0.02 x 10~ 6 and for region 3 even lower at 7.00±0.03 x 
10~ 6 . The C2 component gas in region 3 has a nitrogen abundance 
that is higher than its CI counterpart by a factor of 2. As a result, 
there is a peak in the N/O ratio within UM 448 located at region 2, 
which can be seen in the bottom panel of Fig. 9. For CI, regions 1 
and 3 agree with an average of log(N/0) ~ -1.23±0.04, whereas 
region 2 has an average (N/O) ratio of ~ -1.03±0.05. The map of 
log(N/0) shows that this peak actually reaches ~ —0.80 at the pre- 
cise location of region 2 (as can be seen from the Ha contours). 
The N/O ratio of the (broad) C2 component gas in region 3 (~ - 
1.26±0.13) is in good agreement with the average N/O of CI gas in 
that region. The long-slit study by IT98 measured an N/H ratio of 
9.5 ±2.0 x 10 -6 , in agreement with the averages of regions 1 and 2 
derived from the N/H map shown in Fig. 9. The average CI N/O ra- 
tios for region 3 is slightly higher (by ~0.2 dex) than those of other 
galaxies of similar oxygen metallicity (see Lopez-Sanchez & Es- 
teban 2010, their figure 11). They are, however, consistent within 
the uncertainties. Region 2, however, has a higher value than the 
average value of ~ —1.4 ±0.2 expected for its metallicity, even 
when taking the uncertainties into consideration, for both the av- 
eraged value of region 2 and on a spaxel-by-spaxel basis. For this 
region we do observe the factor of ~ 3 nitrogen excess reported by 
Pustilnik et al. (2004). 

The Ne/H and S/H abundance distribution throughout UM 448 
for the narrow (CI) emission line component show only minimal 
abundance variations in relation to the Ha emission. Ne/H values 
averaged over each region range from ~2.06-6.55 x 10 -5 , with the 
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Figure 9. UM 448 nitrogen, oxygen and log(N/0) abundance maps for velocity components CI and C2, with integrated light Ha emission line contours 
overlaid in dashed line. North is up and east is to the left. 



highest Ne/H abundance aligning with region 2. The Ne/O ratios 
of region 1 and 3 are typical for those expected by BCGs at simi- 
lar metallicities; the averaged region values lie within the reported 
range of log(Ne/0) ~ -0.9 to -0.7 (Izotov et al. 2006). However, 
the Ne/O ratio for region 2 is ~0.5 dex higher than the expected 
Ne/O value. Calculations of Ne/H were also made using the Ne 2+ 
ICF of Izotov et al. (2006), however the Ne/O ratio for region 2 
remained ~0.2 dex higher than the expected value. The S/H distri- 
bution across UM 448 is relatively constant for the CI and C2 gas 
components, with each region displaying averages within the range 
of ~ 4.7 to 7.4 x 10~ 6 . The CI S/O ratio for all regions are slightly 
higher than expected for blue compact galaxies, log(S/0)~ —1.9 



to —1.5 (Izotov et al. 2006), and the S/O value C2 component in 
Region 3 lies at the upper limit of this average. 



5.2 Comparison of spatially resolved and global spectra 
abundances 

Also listed in Table 3 are the T a , N c , ionic and elemental abun- 
dances derived from summed spectra over the entire galaxy (using 
line fluxes listed in Table 2). This particular set of results allows us 
to assess if global galaxy spectra can reliably represent the physi- 
cal properties of its ISM. For example, values derived from long- 
slit observations may contain luminosity-weighted fluxes and con- 
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Figure 10. [O III] A4363 profile observed in the spectra summed over re- 
gion 3, displaying a narrow (CI) plus broad (C2) component Gaussian pro- 
file. Flux, FWHM and velocity parameters for each component are listed in 
Table A3 



sequently erroneous ratios, or large apertures may include a mix- 
ture of gas with different ionization conditions and/or metallicities. 
This question was addressed by Kobulnicky et al. (1999) to assess 
whether global emission-line spectra can be depended upon to reli- 
ably recover the chemical properties of distant star-forming galax- 
ies. Overall, they find that when using direct-method techniques, 
global spectra can give an accurate representation of the oxygen 
abundances derived from individual H II regions, but only after ap- 
plying a correction of A(OfH)— 0.1 to the global spectra results. 
In addition to this, Pilyugin et al. (2012) found that oxygen abun- 
dances were systematically underestimated (and N/O consequently 
over-estimated) when using global SDSS (Sloan Digital Sky Sur- 
vey) spectra of H II regions with different physical properties. 

For UM 448, the temperature derived from the integrated spec- 
trum is in agreement with the flux-weighted mean T c across the 
three regions (14010 ± 2100 K), suggesting that the integrated 
[O III] ratio may have been more heavily weighted by areas of 
stronger [O III] A4363 emission. The oxygen abundance from the 
integrated spectrum, 12+log(O/H)=8.34±0.07, is higher than all 
CI -component values, suggesting that it may be weighted towards 
the high O/H abundance seen in the C2-component of region 2. Be- 
ing +0.18 dex higher than the flux-weighted average of the three 
regions (12+log(O/H)=8.16±0.18), the globally-derived (i.e. de- 
rived from the integrated spectrum) O/H is not fully representa- 
tive of the emission from the whole galaxy, even after applying the 
correction suggested by Kobulnicky et al. (1999). Nitrogen too, dis- 
plays a higher abundance when derived from the summed spectrum 
rather than the regional averages, and results in a (log) nitrogen-to- 
oxygen ratio of —1.15, which is approximately the expected value 
at this metallicity. The Ne/H abundance ratio from the integrated 
spectrum is almost a factor of 2 higher than the flux-weighted aver- 
age of the three regions (~2.14x 10 -5 ), whilst the log(Ne/0) ratio 
is in good agreement. The same applies for the S/H abundance, al- 
though less severe, where the globally-derived S/H ratio is a factor 
~1.5 higher than the flux-weighted average S/H across the three 
regions, but the log(S/0) value is in agreement between the two 
methods. 

Comparisons can also be made with the results from a long- 
slit study by IT98, which utilized summed spectra over a 2" x 7" 
region 'across the brightest region of the galaxy'. From this spec- 
trum they derive an T e of 12100±500 K (in line with the CI com- 
ponent values of our region 1), O, N, Ne and S, abundances that are 
systematically lower than those derived from our summed spec- 
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Figure 11. NTT-SuSI2 combined R, V and S-band images of UM 448 
with contours showing the continuum near Ha (top panel) and integreated 
Ho emission (bottom panel; see Fig. 3 for maps of the separated emis- 
sion line components). North is up and east is to the left. Ha contours are 
shown for the range 35.0 to 128.6 in steps of 22.0, in flux units of 10~ lfi 
ergs s — 1 cm -2 arcsec - 2 . The continuum peak region is represented by 
the white-dashed box; spectra were summed over this region to investigate 
stellar absorption features and averages in EW(H/3). 



trum and abundances relative to oxygen that are in relatively good 
agreement within the uncertainties. However, without knowing the 
exact positioning of this slit, a truly meaningful comparison cannot 
be made. 

Overall it can be said that whilst abundance ratios of heavy el- 
ements can be reasonably well obtained by summed spectra and/or 
long-slit observations, the ionic and elemental abundances relative 
to hydrogen are not well reproduced; this may be because the for- 
mer properties do not have a very severe sensitivity to T c biases, 
while the latter do. This is in addition to the fact that, in contrast 
to integral field spectroscopy, global properties derived from light- 
averaged methods cannot resolve abundance variations and any po- 
tential sites of enrichment are being averaged out. Of course, un- 
certainties will still exist within IFU observations due to averaging 
along the line of sight of individual spaxels; this is something that 
can only be remedied with the development of higher spectral res- 
olution IFUs. 



6 STELLAR POPULATIONS 
6.1 Age of the Stellar Population 

Figure 1 1 shows a combined NTT-SuSI2 R, V and B-band image 
that shows the stellar continuum within UM 448. Overlaid on the 
top panel are the contours of the continuum emission in the Ha 
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8.02± 0.13 


8.42± 0.10 


O. 1 \J.\J 1 


z/z 


0.32± 0.02 


0.22± 0.02 


0.20± 0.03 


0.5 1± 0.05 


0.42 ± 0.03 


N T /« T X 10 


4.19± 0.14 


/ AO 1 f\ r\ Q 

6.08± 0.28 


3.49± 0.21 


8.45± 0.85 


8.74 ± 0.76 


ICF(N) 


2.24 


1.70 


2.00 


1.70 


1.77 


xt /T_r w 1 n6 
JN/ri X 10 


9.40± 0.00 


10.35± 1.40 


7.00± 1.24 


14.40± 2.53 


1 C A A I 1 n 1 

13.44 ± i.yi 


log(N/0) 


-1.24± 0.05 


-1.03±0.09 


-1.18± 0.12 


-1.26± 0.14 


-1.15 ±0.10 


Ne 2 +/_f/+ xlO 5 


2.10± 0.12 


2.70± 0.5 


1.03± 0.13 




3.16 ±0.63 


ICF(Ne) 


1.80 


2.42 


2.00 




2.30 


Ne/H X 10 5 


3.78± 0.31 


6.55± 0.1.46 


2.06± 0.41 




7.27 ± 1.59 


log(Ne/0) 


-0.63± 0.10 


-0.23± 0.25 


-0.71± 0.24 




-0.47 ± 0.20 


S+/H+ XlO 6 


1.11± 0.04 


0.92± 0.05 


1.06± 0.06 


1.43± 0.09 


1.87 ±0.17 


ICF(S) 


1.06 


1.02 


1.04 


1.02 


1.03 


S/H X 10 6 


7.23± 0.45 


4.75± 0.63 


6.30± 1.05 


7.35± 1.06 


9.93 ± 0.92 


log(S/0) 


-1.35± 0.08 


-1.37±0.17 


-1.22± 0.21 


-1.56± 0.18 


-1.34 ±0.12 



Table 3. Ionic and elemental abundances for UM 448, with the corresponding T e and 7V e used in their derivation. Columns 2-5 refer to abundances derived 
from maps, weight-averaged over regions defined in Fig. 2. Column 6 refers to T e , N e and abundances derived from spectra integrated over UM 448 
(corresponding fluxes are listed in Table 2). All ionisation correction factors (ICFs) are taken from Kingsburgh & Barlow (1994). S/H ratios were calculated 
using the ICF(S) values listed along with a predicted value of S 2 + (see text for details). 



spectral vicinity (~6 A windows either side of the line) that are 
in good agreement with the morphology of the broad-band images. 
In the lower panel we show the same broad-band images with the 
Ha integrated-light emission line contours. Significant differences 
between the morphology of the ionised gas emission and that of 
the stellar continuum are clearly apparent, especially in the top, 
western 'arm'. A strong continuum peak can be seen, which is not 
mirrored by the Ha emission, suggesting that a strong stellar com- 
ponent exists within this region. From the broad-band colour im- 
age shown in Figure 1 it can be seen that this peak in continuum is 
red in colour, signifying an older stellar population than that seen 
in the bluer, eastern 'arm'. Such a distinct, two component stellar 
structure, in conjunction with the difference in radial velocity be- 
tween them (see Fig. 3b), supports the interpretation that UM 448 
is a merging system comprised of two galaxies, with very differ- 
ent ages. Spectra summed over this continuum peak region (rep- 
resented by the white-dashed box overlaid in Fig. 1 1 -lower panel) 
do not reveal any stellar absorption features around the Ha, H/3, 
H 7 or HS lines, however Ca II H, K (AA3970, 3934) interstellar 
absorption features were detected. 

The age of a young population in a galaxy can be estimated us- 
ing hydrogen recombination lines, since they provide an estimate of 
the ionizing flux present, assuming a radiation-bounded HII region 
(Schaerer & Vacca 1998). In particular, the equivalent-width (EW) 
of H/3 can be used as an age indicator of the ionizing stellar popu- 
lation at a given metallicity. A map of EW(HP) for UM 448 can be 
seen in Figure 12. Whilst there are no clear peaks in EW(Hf3), the 
main region of decreased EW correlates with the largest peak in 
continuum emission (red contours) and that corresponds to region 
3. 

Following the method outlined by James et al. (2010), we can 



use this map, in conjunction with the metallicity map described 
in Section 5, to estimate the age of the most recent star forming 
episodes throughout UM 448, by comparing the regional observed 
average EW(H/3) values with those predicted by the spectral syn- 
thesis code STARBURST99 (Leitherer et al. 1999). For the mod- 
els we chose a metallicity 0.2 Z Q , the closest representative to the 
weighted average metallicity of UM 448 (~0.25Z Q , derived from 
regional averages), together with assumptions of an instantaneous 
burst with a Salpeter initial mass function (IMF), a total starburst 
mass of 1 X 10 6 M Q , a 100 M stellar upper mass limit (which ap- 
proximates the classical Salpeter (1955) IMF). Models were run for 
Geneva tracks with "high" mass loss rates and Padova tracks with 
thermally pulsing AGB stars included. For each of these evolution- 
ary tracks, two types of model atmosphere were used; firstly the 
Pauldrach-Hillier (PH) and secondly Lejeune-Schmutz (LS). The 
latter was chosen because it incorporates stars with strong winds, 
which would be representative of the WR population within the 
galaxy (if present). The stellar ages predicted by each model and 
the observed average EW(H/3) within each peak emission region 
in UM 448 were derived from Figure 13 and are listed in Table 4. 

The difference between the ages predicted by the Geneva 
and Padova stellar evolutionary tracks is relatively small, with the 
Geneva tracks predicting lower ages by up to 23%. The differ- 
ence in ages predicted by the PH and LS atmospheres are smaller 
still, with LS predicting lower ages by ~5%. We cannot com- 
ment on which atmosphere model is more appropriate given that 
the existence of WR stars within UM 448 is not confirmed by the 
current study, hence we adopt average stellar ages from the four 
model combinations. The current (instantaneous) star formation 
rates (SFRs) based on the Ha luminosities were calculated fol- 
lowing Kennicutt (1998) and are given in Table 4. SFRs corrected 
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EW (H/S) / A 

Figure 12. Map of the equivalent width of H/3 accross UM 448. Overlaid 
red lines are contours of H/3-region continuum; black lines are H/3 flux 
contours from Fig. 2. 



Figure 13. EW(Hf3) as a function of age, as predicted by the STAR- 
BURST ^ code for metallicities of 0.2Zq using a combination of Geneva or 
Padova stellar evolutionary tracks and Lejeune-Schmutz (LS) or Pauldrach- 
Hillier (PH) model atmospheres. The observed average J?W(H/))'s for each 
star forming region in UM 448, along with the peak in continuum high- 
lighted in Fig. 1 1 (lower panel), are overlaid (solid line). 



for the sub-solar metallicities are also given, derived following the 
methods outlined by Bicker & Fritze-v. Alvensleben (2005). 

Regions 1-3 contain ionizing stellar populations with aver- 
age ages of 4.6, 5.5 and 6.0 Myr, respectively, and the region of 
peaked continuum (Fig. 11) has an average age of 8.2 Myrs (Ta- 
ble 4). since Region 1 shows the highest ionizing flux, we would 
expect it to contain the youngest stellar population. It is notable 
that the ages of the stellar population increase through Regions 1-3 
as one moves away from the blue continuum peak (Fig. 1). How- 
ever, the SFR does not show a smooth decrease in the same direc- 
tion, decreasing slightly for region 2 before increasing again for 
region 3. Interestingly, and not surprisingly, the main red-colored 
stellar continuum peak (Figs. 1 and 11) displays an increase in age 
of almost 4 Myrs from the blue continuum peak (region 1) and a 
much decreased SFR compared to the three main regions of star- 
formation. Using narrow-band Ha imaging, Dopita et al. (2002) 
derived a total SFR of 6.6 M©/yr, which is in good agreement with 
that derived from our total Ha flux, 5.0 M0/yr. Accounting for its 
~0.2Zotnetallicity, this translates to 2.8 M©/yr. The evidence sug- 
gests that UM 448 is made up of two separate bodies undergoing a 
merger; the first is sampled by our aperture regions 1 and 2, where 
relatively recent bursts of star-formation have occurred. The second 
body with which it is merging may be an older galaxy (represented 
by the red continuum peak) sampled by aperture region 3, within 
which the collision could have triggered star-formation, hence an 
increased SFR but overall age that is representative of an aggregate 
old+young population. 

6.2 WR stars and N-enrichment 

Contrary to the conclusions of previous long-slit studies (Masegosa 
et al. 1991; Schaerer et al. 1999, and references therein) and in 
agreement with Guseva et al. (2000), our FLAMES IFU spectra 
do not reveal evidence of 'strong' WR emission bands at ~4640- 
4686 A. Unfortunately the C IV 5801,5812 A feature is not included 
in our wavelength coverage, although these features are usually 
significantly weaker than the ~4686 A feature. We find no evi- 
dence of a WR population in our single spaxel spectra and only 
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Figure 14. Sections of continuum-subtracted FLAMES IFU spectra, show- 
ing the blue 'WR bump' region within summed spectra over the entire 
galaxy and separate aperture regions 1, 2 and 3 (see Fig. 2). There are min- 
imal detections of WR emission features within both the spectra summed 
over the entire galaxy and also each of the individual aperture region spec- 
tra. Overlaid are the estimates of 4000±1100, 900±600, 300±600 and 
1300±400 WN (5-6) stars (solid red line), for the entire galaxy and regions 
1—3, respectively, using the templates of Crowther & Hadfield (2006), along 
with the Ict errors (dashed red line). Vertical dashed lines indicate the loca- 
tion of [Fe III] 4658 A and a weak, narrow, He II emission line at 4686 A. 
The dashed horizontal lines in the plots show the zero flux levels. 

minimal evidence in our spectra summed over the entire galaxy. 
When searching for WR features one must consider the width of the 
extraction aperture, which can sometimes be too large and dilute 
weak WR features by the continuum flux (Lopez-Sanchez & Este- 
ban 2008). This has indeed been the case for several other BCGs, 
where previous reports of detections of WR features in long-slit 
spectra were not confirmed by IFU spectra (e.g., James et al. 2010; 
Perez-Montero et al. 201 1). We therefore also examined the spectra 
summed over each individual region (Figure 2), which are shown 
in Fig. 14. 

In order to determine the number of WR stars from the 



14 B. L. James et al. 

Table 4. Age of the latest star formation episode and current star formation rates. 







Starburst aj 


jes (Myr) 




Model 


Reg. 1 


Reg. 2 


Reg. 3 


Cont m Peak 


Padova-AGB PH 


5.06± 0.24 


5.86+ 0.13 


6.36+ 0.48 


8.61+ 0.74 


Padova-AGB LH 


4.96± 0.24 


5.66+ 0.13 


6.01+ 0.33 


8.51± 0.79 


Geneva-High PH 


4.41+ 0.27 


5.31± 0.16 


5.91+ 0.56 


7.91+ 0.74 


Geneva-High LS 


3.86+ 0.48 


5.11± 0.16 


5.61+ 0.50 


7.76+ 0.82 






Star formation rates (Mq yr 




SFR(Ho) a 


1.49±0.03 


0.75+0.02 


0.90±0.02 


0.34±0.01 


SFR(Ha) 6 


0.82±0.02 


0.41±0.01 


0.50+0.01 


0.18±0.01 



a Derived using the relationship between SFR and L(Ha) from Kennicutt (1998). 

b Corrected for sub-solar metallicities following Bicker & Fritze-v. Alvensleben (2005). 



strength of the 4640-4690 A feature for the summed spectra over 
regions 1-3, the local continuum was subtracted. A linear interpo- 
lation to the continuum from both sides of the 4640-4686 A region 
was employed based on the mean flux in a blue region (rest wave- 
lengths 4609 - 4640 A, avoiding the [Fe III] emission line and pos- 
sible C III absorption at 4647 A) and a red region (4711-4742 A). 
Utilizing an LMC WR (WN5-6) spectral template from Crowther 
& Hadfield (2006), the best match of the scaled WR star flux to 
the flux over the observed broad 4686 A region was determined by 
minimizing \ 2 over the peak of the WR feature, typically 4670 - 
4900 A. Since errors were available on the flux points these were 
used to generate a series of Monte Carlo realizations of each re- 
gion spectrum and the WR template flux was best fitted to each. 
The distribution the WR template flux to each realization was fitted 
by a Gaussian enabling the ler errors to be estimated. The flux of 
the WR template was then converted into the number of WR stars 
allowing for the distance of UM 448. We estimate 900 ± 600, 300 
± 600 and 1300 ± 400 WN stars for regions 1-3 respectively. A 
total population of 4000 ± 1 100 WN stars was estimated applying 
the same procedure for the summed spectrum for the entire galaxy. 
These error estimates are based purely on the flux errors and do not 
take into account any systematic error arising from the, necessarily 
simplified, removal of a linear continuum under the WR line.The 
WR template fits and the error ranges for the summed spectra over 
the entire galaxy and the three regions 1-3 are shown in Fig. 14. 
However, considering the large uncertainties, the existence of a WR 
population within UM 448 is still considered to be speculative and 
as a result, we cannot conclude that the Alog(N/0)~ +0.6 dex ex- 
cess observed in region 2 is due to nitrogen enrichment from WR 
stars. 

However, following the methodology outlined in Perez- 
Montero et al. (201 1), we can use the size of the N-enriched region 
in UM 448 (~1.6" x 1.6" = 0.6 kpc 2 ) in region 2 as a scale length 
to assess whether the reported WR stars could be the actual source 
of this N-pollution. Using modeled chemical yields at ~O.2Z , 
those authors study the effects of elements ejected into the ISM via 
winds as a function of time and pollution radii (see Perez-Montero 
et al. 201 1, section 4.3 for a full description of the model parame- 
ters). After ~4 Myr, the appearance of WR star winds increases the 
N/O ratio, where the value of N/O ratio only reaches that of region 
2 at radii much lower (e.g., r = 10 or 100 pc) than those mea- 
sured from our IFU observations. Therefore, although WR stars 
have been detected in previous long-slit studies, and are tentatively 
detected here, the high N/O ratio observed within region 2 is not 
likely produced by the pollution from stellar wind ejecta coming 



from WR stars. Instead, this enrichment may be related to other 
global processes occurring within UM 448, namely the likely in- 
teraction and/or merger between two bodies as inferred from the 
SuSI2 images and radial velocity maps. Studies suggest that in- 
teracting galaxies fall >0.2 dex below the mass-metallicity rela- 
tion of normal galaxies due to tidally induced large-scale inflow of 
metal-poor gas to the galaxies' central regions (e.g. Kewley et al. 
2006; Michel-Dansac et al. 2008; Peeples et al. 2009; Rupke et al. 
2010; Torrey et al. 2012). In addition, models by Koppen & Hensler 
(2005) have shown that a rapid decrease in oxygen abundance dur- 
ing an episode of massive and rapid accretion of metal-poor gas 
is followed by a slower evolution which leads to the closed-box 
relation, thus forming a loop in the N/O-O/H diagram. Such a sce- 
nario, in conjunction with metal-rich gas loss driven by SN winds, 
has been suggested by Amorin et al. (2010) to explain the system- 
atically large N/O ratio found within green pea galaxies. 



7 SUMMARY & CONCLUSIONS 

We have analysed FLAMES-IFU integral field spectroscopy of the 
BCG UM 448 and studied its morphology by creating monochro- 
matic emission line maps. This system shows signs of interaction 
and/or perturbation and it is possibly undergoing a merger. Based 
on the emission line morphology, we defined three regions of high 
Ha surface brightness which we use to derive error-weighted av- 
erages of the physical and chemical conditions across UM 448. 
Region 1 is to the east and corresponds to the tidal tail, region 3 
extends towards to the upper-west and region 2 lies between them. 
UM 448 exhibits complex emission line profiles, with most lines 
consisting of a narrow, central component and an underlying broad 
component, and a third, narrow blue-shifted component is also de- 
tected in Regions 2 and 3. 

Abundance ratio maps of O/H, N/H, Ne/H and S/H were cre- 
ated based on the direct method of estimating electron tempera- 
tures from the [O III] A4363/A4959 line ratio. Minimal variance in 
oxygen abundance are seen across UM 448 for the narrow velocity 
component, with an average abundance of 12+log(O/H)=8.15±0.1. 
However, in an area of UM 448 where the [O III] 4363 A line is re- 
solved into a narrow+broad component, an ~0.4 dex increase in 
oxygen abundance attributable to a decrease in T c of ~5000 K is 
found in the broad line region. We detect the previously reported 
high N/O ratio, with Alog(N/0)~ +0.4, in our aperture region 
2 only, spatially coincident with a peak in nitrogen abundance. 
Abundances derived from an integrated spectrum of UM 448 are 
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not a reliable representation of the ionic and elemental abundances 
throughout UM 448, which may be a consequence of the variance 
in T c throughout the nebular emission line components of UM 448 
detected in these integral field spectra. However, abundances rela- 
tive to oxygen from the two methods (i.e. map-based averages vs. 
integrated spectra) are consistent within the uncertainties. These 
results are significant when assessing the reliability of abundances 
derived from long-slit and integrated spectra of both nearby and 
high-redshift star-forming galaxies that do not contain homoge- 
neous ionization conditions and/or metallicities. 

Continuum images of UM 448 display a rather different mor- 
phology to those of the emission lines, with a peak in stellar light 
situated in between the second and third custom-aperture regions. 
The colour of this continuum peak (from SuSI2 broad-band im- 
ages) shows signs of a red, older stellar population in contrast to 
a younger, bluer stellar population situated in the eastern body of 
UM 448. Star formation rates were computed from the Ha line lu- 
minosities within each sampled star forming region and the ages of 
the last major star formation episode were estimated by fitting the 
observed Balmer line equivalent widths with STARBURST99 mod- 
els. The two merging sections of UM 448 have SFRs that differ by a 
factor of ~5 and starburst episodes separated by 4 Myr, providing 
evidence of a distinct age gradient between the two components. 
The latest major star formation event took place ~4.5 Myr ago in 
the largest merging component, and has been producing stars at a 
rate of 0.8 M yr" 1 . 

Finally, we do not detect the previously reported 'strong' WR 
signatures in UM 448, on a single-spaxel, regional or global scale 
(Masegosa et al. 1991; Schaerer et al. 1999). Using the 4686 A fea- 
ture we tentatively estimate that ~4000 WN stars may be present 
within UM 448. In consideration of the physical size of increased 
N/O seen within UM 448 (~0.6 kpc 2 ) we find that WR stars 
alone cannot be responsible for the region of N-enrichment within 
UM 448. Thus being, another chemical evolution scenario (e.g. in- 
flow of metal-poor gas due to interaction or merger scenario, as the 
morphology and dynamics of UM 448 suggest) is required to ex- 
plain this enhancement. Since there is no strong evidence in support 
of these scenarios (although broad component emission is sugges- 
tive of an outflow, the N-enrichment is seen in the narrow com- 
ponent only), more observations, at higher sensitivity, are needed 
of this intriguing system. In relation to other high N/O BCGs, this 
chemodynamical study supports its membership within the anoma- 
lous subset and, along with other spatially resolved spectroscopic 
studies of high N/O galaxies, demonstrates the benefits that this 
method can have in both isolating the enrichment and aiding our 
understanding of the processes that cause it. 
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APPENDIX A: REGIONAL OBSERVED FLUXES AND 
LINE INTENSITIES 

Tables A1-A3 list the observed fluxes and de-reddened line intensi- 
ties derived from summed spectra over the regions defined in Fig. 2. 
Fluxes are given for each separate velocity component (C1-C3), 
along with its de-reddened line intensity and FWHM. 

This paper has been typeset from a TgX/ ETgX file prepared by the 
author. 
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CI 






C2 






Fa 


l \ 


FHWM 


Fa 


Ia 


FHWM 


[On] A3727 


93.39±1.84 


108.46±3.62 


77.07±0.48 


69.57±2.33 


100.49±9.06 


167.96±1.65 


[0 11] A3729 


110.36±2.38 


128.17±4.43 


77.02±0.48 


129.04±3.49 


186.39± 16.39 


167.81±1.64 


[Ne 111] A3868 


8.69±0.35 


9.94±0.48 


58.50± 14.57 


35.36±3.74 


49.28±6.60 


169.82± 15.62 


H8+He I A3888 


10.69±6.14 


12.22±7.02 


70.76±17.14 


25.16±8.63 


34.91 ±12.31 


172.18±38.13 


H<5 


18.53±1.38 


20.61±1.63 


80.58±2.20 


32.73±2.02 


42.51 ±4.25 


145.87±2.70 




A'i 7S_|_T 1 A 


AH 1 ZA 


Sfi 1^-1-1 


Aft 1 7_1_ 7 77 
■+U. 1 / ztZ. i L 




1 7*> ^7 

1 / J.Jji'+.jZ 


[0 ill] A4363 


2.93±0.27 


3.15±0.30 


104.09± 10.38 








He l A4471 


3.40± 1.01 


3.60±1.07 


74.58±14.14 


5.02±1.27 


5.76±1.52 


236.59±52.09 


H/3 


100.00±2.23 


100.00±3.09 


81.32±0.69 


100.00±3.01 


100.00±7.31 


189.73±2.19 


[0 ill] A4959 


91.61±2.46 


90.34±3.08 


76.61±0.82 


100.14±3.49 


96.75±7.14 


182.57±2.65 


[Nil] A6548 


15.15±0.57 


12.58±0.51 


69.33±1.26 


14.51±0.76 


9.19±0.64 


156.65±4.10 


Ha 


349.49± 18.66 


291.23±16.13 


104.33±0.35 


467.65± 10.91 


298.72± 15.36 


235.81 ±1.00 


[N 11] A6584 


56.17±1.15 


46.50±1.17 


77.19±0.65 


37.07±1.15 


23.31±1.28 


194.92±3.40 


He l A6678 


0.50±0.16 


0.41±0.14 


30.78±6.20 


9.16±0.29 


5.65±0.31 


87.88±1.67 


[S 11] A6716 


37.41±0.82 


30.63±0.80 


71.74±0.55 


40.41±1.11 


24.73±1.28 


155.42±1.32 


[S 11] A6731 


29.24±0.63 


23.95±0.62 


71.59±0.55 


0.91 ±0.02 


0.56±0.03 


155.08±1.32 


c(H/3) 




0.25 ± 0.01 






0.62 ± 0.03 




F(H/3) 




5.30 ± 0.08 






3.83 ± 0.08 





Table Al. UM 448 regional fluxes and de-reddenened line intensities (both relative to F(H/3)=I(H/3)=1()0) and FWHMs (in km s _1 ) for summed spectra over 
region 1 (as defined in Fig. 2). Line fluxes were extinction-corrected using the c(H/3) values shown at the bottom of the table, calculated from the relative Ha, 
H/3 and H7 fluxes. F(H/3) is in units of X 10~ 14 ergs s^ 1 cm~ 2 arcsec~ 2 . 







CI 






C2 






Fa 


Ia 


FHWM 


Fa 


Ia 


FHWM 


[Oil] A3727 


135.44±4.93 


165.56± 11.27 


88.05±1.05 


73.89±3.39 


97.79±7.94 


197.23±3.84 


[Oil] A3729 


160.52±6.28 


196.21±13.65 


87.98±1.05 


117.41±4.18 


155.40±11.78 


197.06±3.84 


[Nem] A3868 


30.52±5.22 


36.59±6.59 


102.24± 10.23 


14.84±3.93 


19.11±5.21 


323.00±92.41 


H8+He 1 A3888 


8.18±7.09 


9.78±8.49 


60.54±25.61 


25.41±5.01 


32.61 ±6.77 


152.59±21.87 


H5 


14.24±3.86 


16.43±4.54 


80.17±7.25 


27.99±2.90 


34.16±4.14 


139.71±3.95 


H 7 


42.40±3.00 


46.84±4.10 


94.70±4.38 


41.56±2.73 


47.76±4.25 


170.55±7.05 


[O III] A4363 


3.10±0.54 


3.41±0.62 


124.72±24.87 








Hel A4471 


1.43±1.89 


1.54±2.04 


56.25±34.72 


4.60±1.30 


5.11±1.47 


140.16±29.63 


H/3 


10O.O0±3.78 


100.00±5.94 


91.60±1.35 


100.00±3.80 


100.00±6.55 


200.23±3.45 


[Om] A4959 


102.17±3.89 


100.26±5.89 


86.57±1.55 


67.05±3.16 


65.31 ±4.59 


191.52±4.09 


[Nil] A6548 


23.57±1.59 


18.37±1.36 


81.05±2.43 


12.88±1.07 


9.09±0.82 


165.02±5.90 


Ha 


370.44± 11.04 


290.03± 12.58 


103.56±0.57 


415.93± 12.45 


295.55± 13.99 


278.17±2.12 


[N 11] A6584 


97.70±2.92 


75.84±3.27 


92.68±0.88 


31.35±1.13 


22.01±U2 


253.20±4.52 


He 1 A6678 


9.41±0.33 


7.23±0.34 


102.59±2.28 








[S 11] A6716 


53.03±2.02 


40.55±1.97 


85.70±1.06 


29.47±1.24 


20.26±1.11 


180.53±2.77 


[S 11] A6731 


43.70±1.59 


33.42±1.58 


85.52±1.06 


18.79±0.91 


12.92±0.77 


180.14±2.76 


c(H/3) 




0.34 ± 0.02 






0.48 ± 0.02 




F(H/3) 




1.95 ± 0.05 






2.72 ± 0.07 





Table A2. Same as for Table Al, for summed spectra over region 2 (as defined in Fig. 2). 
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CI 






C2 






Fa 


Ia 


FHWM 


Fa 


l \ 


FHWM 


[0 II] A3727 


37620±7.41 


506.95±52.62 


199.37±0.73 








[0 II] A3729 


487.54±9.55 


656.99±68.17 


212.46±15.14 








[Nem] A3868 


65.57±4.66 


85.82±10.52 


213.38±19.76 








H8+He I A3888 


55 24±5.09 


72.05±9.78 


221.94±2.30 








HS 


9.51 ±2. 11 


11.76±2.85 


95.75±9.88 


26.11±1.09 


32.67±1.87 


157.59±5.69 




.1 1 . o ) It J . J J 


Af. 9^-1-^ 7T 
H-U.ZJ It J. / J 


J 1 D.L JZt — O.Uo 






TJfi lf,-\-00 96 
ZZU. /O^ — l. — U 


[0 III] A4363 


1.45±0.84 


1.67±0.97 


74.34±37.30 


5.33±0.07 


6.19±0.24 


455.65±67.72 


He I A4471 


13.64±1.33 


15.25±2.01 


226.35±22.75 








H/3 


100.00±2.71 


100.00±8.56 


118.95±8.93 


100.00±0.96 


100.00±3.46 


233.89±7.09 


[0 III] A4959 


97.79±2.17 


95.09±7.82 


94.79±3.52 


89.70±0.97 


87.09±2.98 


255.45±4.99 


[N II] A6548 


60.50±2.56 


41.76±2.91 


188.97±2.54 








Ha 


414.10±8.58 


287.88±17.14 


111.69±4.90 


417.95±3.94 


284.58±7.03 


249.28±5.71 


[N II] A6584 


132.16±2.76 


90.71 ±5.33 


162.09±7.08 


35.66±3.51 


23.95±2.42 


333.30±17.45 


He I A6678 


14.28±0.71 


9.66±0.71 


200.80±7.71 








[S II] A6716 


176.68±4.77 


1 18.61 ±7. 10 


193.47±0.69 








[S II] A6731 


135.17±4.06 


90.75±5.56 


193.01 ±0.69 








c(H/3) 




0.51 ± 0.04 






0.53 ± 0.01 




F(H/3) 




1.67 ±0.03 






3.54 ± 0.02 





C3 




Fa 


Ia 


FHWM 


[On] A3727 








[O II] A3729 








[Nem] A3868 








H8+Hel A3888 








HS 








H 7 


41.93±3.60 


47.58±5.45 


106.30±20.01 


[O III] A4363 








He I A4471 








H/3 


100.00±5.80 


100.00±8.88 


84.49±6.54 


[O III] A4959 


129.25±8.12 


126.20± 11.47 


86.50±4.99 


[N II] A6548 








Ha 


401.60±16.70 


294.40± 18.32 


86.10±4.71 


[N II] A6584 


77.15±3.50 


55.94±3.60 


91.24±6.78 


Hel A6678 








[S II] A6716 








[S II] A6731 








cW) 




0.43 ± 0.03 




F(H/3) 




0.54 ± 0.02 





Table A3. Same as for Table Al, for summed spectra over region 3 (as defined in Fig. 2). 



